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Abstract 

This paper presents a formulation of lattice fermions applicable 
to ali quark masses, large and small. We incorporate interactions 
from previous light-fermion and heavy-fermion methods, and thus 
ensure a smooth connection to these limiting cases. The couplings in 
improved actions are evaluated for arbitrary fermion mass mg, with- 
out expansions around small- or large-mass limits. We treat both 
the action and external currents. By interpreting on-shell improve- 
ment criteria through the lattice theory's Hamiltonian, one finds that 
cutoff artifacts factorize into the form bn{mqa)[paY'% where p is a 
momentum characteristic of the system under study, s„ is related 
to the dimension of the nth interaction, and hn{mqa) is a bounded 
function, numerically always 0(1) or less. In heavy-quark systems 
p is typically rather smaller than the fermion mass mq. Therefore, 
artifacts of order {mqaY do not arise, even when m^o ^ 1. An impor- 
tant by-product of our analysis is an interpretation of the Wilson and 
Sheikholeslami-Wohlert actions applied to nonrelativistic fermions. 



1 Introduction 



The most promising avenue for a quantitative understanding of nonperturbative 
quantum chr omo dynamics — and other field theories — is via numerical (Monte 
Carlo) integration of functional integrals defined on a lattice [ Q]. Like any nu- 
merical technique this method has uncertainties that must be understood and 
controlied before the resulta are useful. In particular, although the continuum 
theory is defined by the limit of a sequence of lattice theories, the numerical 
calculations are never carried out at the limit. Because the Monte Carlo intro- 
duces statistical errors, the extrapolation to the continuum limit is imperfect. 
The results for physical quantities are consequently contaminated by lattice arti- 
facts. For a practical result, this uncertainty must be smaller than, say, relevant 
experimental uncertainties. 

The way to reduce lattice artifacts is based on the renormalization group [ 

One starts with a general action 

S = ^CnSn, (1.1) 

n 

where the Sn include ali interactions with the desired field content and the ap- 
propriate symmetries. One approach to the continuum limit, which might be 
called brute force, is to choose the c„ in any way that drives the lattice spacing 
to zero. An ideal approach would be to choose the c„ to lie on a renormal- 
ized trajectory [ ^, where there are no lattice artifacts even though the lattice 
spacing a ^ 0. In the space of ali possible actions specified by eq. (|l.l| ), the 
renormalized trajectories lie in a subspace, whose dimension equals the num- 
ber of relevant parameters. Once the relevant parameters have been fixed by 
physics, they and the renormalization scheme determine ali the c„. 

Unfortunately, a renormalized trajectory is mostly of abstract value, because 
on one infinitely many c„ are nonzero. Ali practical schemes, such as blocking 
fields [ ^] or Symanzik improvement [ Q| use criteria such as locality [ ^ or the 
scaling dimension [ ^] to truncate the space of actions. (For an asymptotically 
free theory, such as QCD, these two criteria are not very different.) Furthermore, 
the calculations of the Cn are, in practice, only approximate. For these reasons 
an improved action is only partially renormalized. Nevertheless, any practical 
action can be written 

S = Skt- 6S, (1.2) 

where 5rt denotes (an action on) the renormalized trajectory. Usually the 
truncations and/or approximations used to generate S will also yield estimates 
for the remaining cutoff effects 6S. 

This paper treats massive fermions coupled to non-Abelian gauge fields. The 
relevant couplings are the fermion masses and the gauge coupling. So the renor- 
malized trajectory takes the form 

5RT(mg/AQCD, AgcDo) = Xlcn("ioa,5o)'S'n, (1-3) 

n 
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where niq denotes the fermion mass,[| Aqcd is the scale characteristic of the 
gauge theory, and the argument AqcDa labels the renormahzation point. The 5„ 
are gauge-invariant combinations of four-component fermion and anti-fermion 
fields (ijj and ■0) and the lattice gauge field (Ufj_). For later calculational conve- 
nience we choose the bare, rather than some physical, fermion mass rn-oa and 
gauge coupling to parameterize the couplings c„. 

As a — > the fermion mass niqU is formally smaller than gQ. (By asymptotic 
freedom g^ = [/Jq log(AQQj-,a^)]^^ as a ^ 0.) It is therefore tempting to expand 
the couplings Cn{mQa, g^) in moa, as in previous analyses [|5[ |6|, 0]. But there may 
be fermions satisfying mg/AqcD ^ 1; the charm, bottom, and top quarks are 
examples in nature. If, in practice, m^a is not small, perturbation theory in m^a 
need not be useful, even though perturbation theory in g^ might be. Indeed, 
this regime includes the charm and bottom quarks at currently accessible lattice 
spacings. 



The static [ ^, |^ and nonrelativistic [ 10, 11, ^] effective theories address 
the problems of heavy fermions. Their restriction to ruq ^ Aqcd implies that 
couplings of interactions between particle and antiparticle states may be chosen 
to vanish, and the remaining interactions in eq. (LI) are organized according to 
a p/ruq expansion. But for some niq ~ 2Aqcd the expansion is no longer useful. 
Furthermore, radiative corrections induce power-law terms, e.g. gQ/{mqa), which 
must be canceled by adjusting the Cn- These terms, which diverge as a — > 0, are 
a reminder that the effective theories are to be used at scales below (large) niq. 
Their presence implies that cutoff effects in the effective theories should be 
removed not by brute force, but by keeping a ~ and constructing actions 



systematically closer to the renormalized trajectory [ 11, 12 



This paper presents a way to encompass both the small and large mass 



formulations. The doubling problem is handled by Wilson's method [ 13|. In 
addition to treating the rnoa dependence of the couplings exactly, the centrai 
idea is to enlarge the class of interactions considered in eq. ( |l.3D to include 
those from both the small moo and the large m-rj/AgcD limits. In particular, 
we do not impose a symmetry between couplings of interactions that would 
be related by interchanging the time axis with a spatial axis. For any m^a, 
the four-momentum of a quark in most interesting physical systems satisfies 
[E{p) — E{0)]a <C 1 whenever pa <C 1. But when niqa > 1 the characteristic 
four-momentum of the physics does not respect time-space axis interchange. 
Under such circumstances it is inconvenient and unnecessary to choose an axis- 
interchange symmetric action]^ 

^We use niq for a quark mass defined by a physical condition and moa for the coupUng 
appearing in the action. 

^Rehnquishing axis-interchange symmetry is common in treatments of heavy quarks with 
momentum-space and dimensionai regulators. It is possible to derive deviations from heavy- 
quark symmetry from the Dirac action [ |l4[ |l^ , while maintaining time-space axis interchange 
invariance as a coroUary to Euclidean invariance. But usually the derivations are easier with a 
nonrelativistic action [ |lo[ |l^] , the so-called the heavy-quark effective theory [ ^ . 
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In the appropriate limits, our formulation of lattice fermions shares proper- 
ties of previous ones. On one hand, at dimension five or less, couplings related 
by axis interchange become identical in the limit rriQa 0: the Wilson action 
and the improved action of ref. [ ^ are recovered. But when moa the mass- 
dependent renormalization leaves lattice artifacts that are proportional to pa 
and AqcDa, not m^a. At higher dimension, however, we retain Wilson's time 



derivative and incorporate "spatial-only" interactions into eq. (1.1). 

On the other hand, for niq ^ Aqcd one can interpret the lattice theory in 
a nonrelativistic light. Indeed, ali members of our class of actions approach a 
universal static limit as rriQa — > oo. For moa large but finite, the 1/m^ corrections 
to the static limit can be recovered systematically, provided the fermion mass is 
defined through the kinetic energy, and provided the general action, eq. (|1.3D , is 
truncated only at dimension 4 + s (or higher) . Unlike previous implementations 
of nonrelativistic fermions, however, our approach crosses smoothly over into the 
regime of tiny lattice spacings, where moa <C 1 even for a heavy quark. Thus, 
after several Cn have been tuned dose to a renormalized trajectory, thereby 
removing the worst lattice artifacts, a little brute force can remove the rest. 

Because we make no assumptions about the ratio of fermion masses to other 
scales, our formulation is especially well suited to fermions too heavy for small 
moa methods yet too light for large mg/AqcD methods. With the actions given 
below one can test whether a given fermion is heavy enough to be treated nonrel- 
ativistically, without resorting to brute- force simulations. A practical example 
might be the charm quark, which has a mass only a few times Aqcdj yet even 
on the finest lattices available today mcho is largish, at least i. 

For a concrete determination of the c„, one must choose a renormalization 



group, a criterion for truncating the sum in eq. (1.3), and a strategy for deter- 
mining the c„ . For illustration we adopt here a Symanzik-like procedure [ Q , 
organizing the interactions by dimension. Carrying out this program to arbitrar- 
ily high dimension would produce a renormalized trajectory of a renormalization 
group generated by infinitesimal changes in a. For simplicity, however, most of 
this paper treats interactions only up to dimension five. Although a nonpertur- 
bative determination of the couplings is possible in principle, this paper makes 
the further approximation of perturbation theory in Qq, that is 

oo 

c„ (moa, c/o) = 9o c^n i'^'Oa) ■ (1-4) 

1=0 

Except for sect. ^ we work to tree level, so we often abbreviate Cn\mQa) by 
Cn("ioa). (Explicit one-loop calculations are in progress [ |l^.) Within these 
approximations we determine the c„ by insisting that on-shell quantities take 



their desired values, as first suggested by Lùscher and Weisz [ 19]. 

One calculational procedure is to work out n-point on-shell Green functions 
via Feynman diagrams and tune them to the continuum limit, to the appropri- 
ate order in p. An example is in sect. 0. Because this strategy is limited to a 
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finite number of quantities, it is nontrivial to assume that other quantities are 
improved too [ ^. An alternative is developed in sect. |5|. Starting from the 
transfer matrix we derive an expression for the fermion Hamihonian, valid (at 
tree level) for states with pa <C 1. Because the Hamiltonian is an operator, im- 
proving it to some accuracy guarantees the improvement to the same accuracy 
of infinitely many c-numbers. We show that by adjusting the couphngs Cn{moa) 
correctly, and ahowing physically unobservable redefinitions of the fermion field, 
one can tune the Hamiltonian to the continuum limit, i.e. to the Dirac Hamil- 
tonian, to the appropriate order in p. 

In addition to equations for the c„, the analysis of the Hamiltonian yields two 
important results. One is that the Wilson time derivative needs no improvement. 
The other is a canonically normalized fermion field that, to the accuracy of 
the improvement, obeys Dirac's (continuum) equation of motion. This field is 
a potent ingredient in calculations of matrix elements involving heavy-quark 
systems; see sect. @. 

Owing to the approximations introduced — the truncation of interactions and 
perturbation theory — some cutoff effects remain. If these errors are small, they 
may be estimated by insertions of 5S in correlation functions. If the series for 
c„(moa, ^o) has been developed to L„-th order 

oo 

^S = Y. E 9oC^^inioa)Sn. (1.5) 

A typical term in 5S distorts an on-shell correlation function by an amount of 
order g'^'c^ (moa)(pa)^", where Sn = dimSn — 4, and l > Ln + 1. (If 5„ is 
omitted from the action altogether, then here l = 0.) The analysis presented 
below shows that the tree-level, lower-dimension Cn\mQa) are well-behaved for 
ali masses. We also show that loop diagrams have the same or more benign 
behavior at large mass. In particular, as moa — > oo the Cn {rriQa) either approach 
a Constant or fall as {mga)~^, for some s < s„. These results provide evidence 
that the higher-dimension Cn{mQa) are well-behaved too. 

In Monte Carlo programs it is customary to parameterize the action by the 
hopping parameter instead of the mass. In this notation the dimension-three 
and -four interactions are written 



n n 

sYl [^n{rs - Ji)Un,i'>Pn+ì + V'n+ì('^s + HWl^i^n 



n 



(1.6) 



Some terms of dimension five — to solve the doubling problem [ ^] — are included 
here too. The relation between the bare mass m^a and the hopping parameters 
Ut and Kg is given below in sect. ^. Eq. ( |1.6| ) illustrates how our program 
subsumes properties of the familiar small-mass and large-mass formulations. 
Imposing axis-interchange invariance would set Ks = nt and = 1, and then Sq 
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reduces to the Wilson action [[T^. Rewriting r^K^ = Cg and setting Kg to zero 
with Cg produces the simplest nonrelativistic action [ ^] . 

This pattern continues for dimension-five interactions. Aside from the Wil- 
son terms in Sq, the other dimension-five interactions are the chromomagnetic 
interaction 

Sb = ^CBKg ^ eijki^nCTijBn-klpn, (1-7) 
n;i,j,k 

and the chromoelectric interaction 

Se = iCEKg ^ IpnmEniilpn, (1-8) 
n;i 

where B and E are suitable functions of the lattice gauge field U, as in sect. |2|. 
The light-quark formalism of refs. [ 0, R] considera the special case cb = ce, 



whereas the heavy-quark formalism of refs. [ 11, 12| sets ce = 0. 

The couplings Vg, ( = Ks/kì, cb, and ce are specific examples of the cou- 
plings Cn in eq. (|1.3|). On the renormalized trajectory they are, therefore, ali 



functions of moa. Sect. M shows how to adjust Vg and C so that the relativis- 



tic energy-momentum relation E = yrn"^ + + (5i?iat is obtained for ali mga. 
With the correct choice, for which C 7^ 1 when rrioa ^ 0, the (tree-level) lattice 
artifact (J-Biat is proportional to p^a^ for m^a <C 1 and p'^a/m^ for m^a » 1. 
Similarly, results in sect. ^ include functions CB{rnQa) and C£;(moa) that reduce 
lattice artifacts in the quark-gluon vertex functions to O(p^a^) for m^a <^ 1, or 
(yet smaller) 0{p'^a/mq) for rriqa » 1. 

In their on-shell improvement program refs. [ |l^, ^ introduced changes of 
variables, or isospectral transformations, to expose redundant interactions. Since 
the coupling of a redundant interaction does not influence physical quantities, 
one can choose it according to theoretical or computational criteria distinct from 
improvement. Sect. ^ examines the isospectral transformations when time-space 
axis-interchange symmetry is not imposed. In our formulation many isospectral 
transformations are exploited to keep the time discretization, and hence the 
transfer matrix, as simple as possible. 

The remaining redundant directions can be classified in the Hamiltonian 
approach developed in sect. In a Euclidean version of the standard Dirac- 
matrix basis, given in sect. ^, matrices are either block diagonal or block off- 
diagonal. Block-diagonal transformations are absorbed into a generalized field 
normalization. On the other hand, block- OjQ^-diagonal transformations (called 
Foldy-Wouthuysen-Tani transformations [ in atomic physics) generate leeway 
in choosing the mass dependence of associated couplings.^ For example, in the 
action So + Sb + Se one may freely choose r^, as long as the choice circumvents 
the doubling problem. 

■^The off-diagonal interactions are precisely the ones usually omitted from nonrelativistic 
formulations, yet their presence in our formulation permits a smooth transition from farge to 
small mass. 
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Our approach breaks down, just as any lattice theory does, when pa is too 
large. Fortunately, the typical momenta and mass splittings of hadronic systems 
usually are bounded; the energy scale around the fermion mass is dynamically 
unimportant. In quarkonia the typical energy- moment um scales are niqV and 
mgv'^/2, i.e. 200-800 MeV for charmonium and 200-1400 MeV for bottomonium. 
Similarly, in light and heavy-light hadrons the typical momentum scale is Aqcd , 
i.e. 100-300 MeV. In some processes, such as a decaying heavy-quark system that 
transfers ali its energy into light hadrons, a large three-momentum \p\ ~ rug does 
arise. Then our formulation and its predecessors ali require further extensions. 
One should appreciate, however, that the breakdown arises not from the large 
fermion mass per se, but from the large momentum of the decay products. 

A by-product of our formalism applies to existing numerical calculations, 
done with axis-interchange invariant actions. For niga^l (and, furthermore, 
mg /Aqcd ^ 1) we derive in sect. |9| a nonrelativistic interpretation of such ac- 
tions. One then sees that, with a proper definition of the fermion mass, any 
action described by Sq + Sb + Se, including the Wilson and Sheikholeslami- 
Wohlert fermion actions, has the lattice-spacing and/or relativistic inaccuracies 
of a typical nonrelativistic action. A practical bonus of the nonrelativistic regime 
is that it is no longer necessary to adjust Hs differently from Kf. In heavy-light 
systems, one may also set = c^;. 

This paper is organized as follows: Sect. ^ introduces some notation, includ- 
ing a form of the action better suited to perturbation theory, and the Dirac- 
matrix convention used in later sections. The isospectral transformation of ref. [ 
^] is reviewed and generalized in sect. ^, to determine which couplings are redun- 
dant. Then, to derive improvement conditions, Feynman-diagram methods are 
discussed in sect. ^, and the Hamiltonian method in sect. ^. With a Hamiltonian 
description of the lattice theory in hand, sect. [^estimates cutoff effects in various 
hadronic systems. Sect. ^ considers perturbations from the electroweak interac- 
tions, needed for the phenomenology of the Standard Model [ 21|. Some of the 
issues beyond tree level are outlined in sect. |^. The relationship of our work, 
for rrig/AgcD ^ 1, to nonrelativistic QCD is pursued in sect. |9|. We discuss a 
few phenomenologically relevant applications more thoroughly in sect. Fi- 
nally, sect. 11 contains selected concluding remarks, and the appendices contain 
various technical details. 



2 Notation 

We shall cali the form of the action in eq. (|l.6| ) the "hopping-parameter form." 
For studying the continuum limit and developing perturbation theory it is useful 
to present a different form. Let us introduce some notation. The lattice spacing 
is a and the site labels are n = x/a. Rescale the fields: 

3/2 

V^n = ^7=Hx) (2.1) 
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and similarly for ipn- The bare mass is 

moa = -^-[l + r,C(d-l)], (2.2) 

2Kt 

where d (= 4) is the spacetime dimension, and (" = Kg/i^t- With these substitu- 
tions the action reads 



5o = mo / i>{x)i,{x) + J [^{xmi + 7o)Doi^ix) - V(x)i(l - 7o)D+ÌJÌx) 
+C / Ì>{x)-r-DTp{x) - iar^C / i>{x)A^^'^tlj{x), 



(2.3) 

where the integrai sign abbreviates J2x • The covariant difference operators 
are conveniently defined via covariant translation operators 

T±fM'>P{x) = U±f,{x)ip{x ± afj), ■4'{x)T±^ = ^{x =F afj)U±^{x =F afi), (2.4) 

where U-^{x) = U^{x — a fi). Then 

L'+V = a-i(ro-i)v, 

Z?o"V' = a"'(l-r_o)V', 

Dixl^ = {2a)-\Ti-T_i)i^, (2-5) 

A(3)V; = a-2^(r, + r_,-2)V', 

define various covariant difference operators and the three-dimensional discrete 
Laplacian. We shall cali the form of the action in eq. (|2.3| ) the "mass forni." 



The temperai kinetic energy in eq. (2.3) is written in a way that does not 
make the temporal Wilson term explicit. Eq. ( |2.3| ) is more convenient, how- 
ever, for constructing the transfer matrix, and for comparing with nonrelativis- 
tic QCD. The spacelike Wilson term, the one proportional to r^, is needed to 
prevent doubling. A convenient choice in computer programs is = 1, but we 
keep it arbitrary, because other choices may have other advantages. 

For constructing the transfer matrix and for examining the nonrelativistic 
limit, a useful representation of the Euclidean gamma matrices is 

satisfying {7^,71/} = 25^^. Another convention that we use is = §[7^,7;^] so 
that (T^j^ = +cj^u- Using eq. (|!|), cjoì = iai and Uij = -Eijk^k, where 
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The following split into upper and lower two-component spinors 



ti; 



X 



(2. 



follows from eq. ( p.6D . This convention is chosen so that (the operators corre- 
sponding to) (/> and x annihilate particle and anti-particle states, respectively. 
With these two-component fields the mass form of the action is 



So = mo J y)^x)(l){x) + xHx)x{x 

+ J [(l)^{x)D^(l){x)+xKx)DoX{x 
-iar^C J [0l'(x)A(3)0(x) + x1'(x)A(3);^(x) 
+ C / \(l)^x)(T-Dx*{x) - x^{x)(T-D(j){x 



(2.9) 



This form of the action exhibits exphcitly that particles and anti-particles are 
treated on the same footing. (Anti-particles transform under the complex- 
conjugate representation of the gauge group, however, so U* appears instead 
of U in the rules (eqs. ( |2.4| )) for covariant translations acting on x.) 

Writing Bn = a'^B{x) and En = o?E{x), the four- and two-component mass 
forms of the chromomagnetic and chromoelectric interactions are 



Sb = -^acBC / Ì'ix)'S-B{x)^{x) 



\acBC J<j)\x)a-B{x)(j){x) - x\x)a* ■B{x)x{x), 



(2.10) 



and 



Se = —^uceC / '^l^{x)o^■E{x)^l^{x) 



\acEC j (l}^{x)(T-E{x)x*{x) + x^{x)<j-E{x)(l){x), 



(2.11) 



respectively. Except in a technical step in sect. ^ we take the "four-leaf dover" 
lattice approximant to the field strength 



1 



F^u{x) = ^ sign(/xi/) Tp,TpT_p_T_p - h.c, 



(2.12) 



fj.=±ii 

u=±u 



as introduced in ref. [ In eqs. ( p.lO ) and ( 2.11| ), Bi = ^eij^Fj^ and Ei = Fqì. 
As defined here, F^i,, Bi, and Ei are anti-Hermitian; similarly we take anti- 
Hermitian gauge-group generators t"^ = —t"". 
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Table 1: Interactions that could appear in the action, with and without axis- 
interchange symmetry (a.i.s.). 



dim 
3 
4 
5 



w/ a.i.s. 



w/o a.i.s. 



3 Redundant Couplings 

Before trying to determine the mass dependence of the couphngs C,, r^, cb, 
and Ce, one should estabUsh which combinations are physical. The fields in 
functional integrals are integration variables, and a change of variables cannot 
change the integrals. Interactions that are induced by changes of variables are 
redundant; their couplings can be chosen with some leeway, dictated by calcula- 
tional or technical convenience, rather than by physical criteria. 

A subtle example of a redundancy in the space of interactions is wavefunc- 
tion (re)normalization, which multiplies the field by a Constant. For fermions, 
for example, it is sometimes convenient to fix the kinetic energy ■^l)p^l; to have 
coefìicient unity, which is the mass form of the action, and sometimes to fix the 
locai term ipip to have coefficient unity, which is the hopping-parameter form. 
But neither interaction is redundant, even though the normalization convention 
drops out of physical quantities. 

Otherwise redundant directions are exposed by redefinitions of the field. In 
the analysis of ref . [ ^] , with axis-interchange symmetry, one considers the trans- 
formation 

where a is chosen so that a{p + m) is "small." After carrying out the transfor- 
mation on the target action f ip{p + m)ijj, one expands the transformed action 
to 0(a). One finds changes in the normalization of the lower-dimension terms 
and the additional interaction a{e+e)'ìpp'^'tp: from the two independent transfor- 
mation parameters, only one combination survives. Hence, of the dimension-five 
interactions listed in Table 0, one (i.e. ipp'^ìp) is redundant, and the other is not. 

On the lattice the nearest-neighbor discretization of p suffers from the dou- 
bling problem. Wilson's prescription adds a nearest-neighbor discretization 
of to eliminate the unwanted states. By the preceding analysis [ us- 
mg instead p^ = D"^ 

~ i^^fiuF^u would not change the spectrum at 0(a). When 
the discretization is chosen to solve the doubling problem, however, the 
interaction does change the spectrum of high-momentum states. Since they 
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communicate with the low-momentum states through virtual processes, lattice 
artifacts proportional to remain. They can be eliminateci with the other 
dimension-five interaction, iipa^uF^^ip, with a coupling proportional to g^. 

Thus, with axis-interchange symmetry there are four interactions up to di- 
mension five, one of which goes with wavefunction normalization (e.g. ippìl^). 
One coupling is redundant, and it can be chosen to solve species doubling 
The other couplings are fixed by the fermion mass (V'V') a physical 
improvement condition {i'ìpCTfj,i,F^^ijj). 

When axis-interchange symmetry is given up, the transformation in eq. ( p.!]) 
should be generalized to 

tph^exp (ea{p + m) + óa-y-Djtl), 

- - \ - \ (3-2) 

■ip h->- tp ex.p lea{p + m) + 6af-Dj. 

Applying this transformation to the target action induces the dimension-five 
interactions listed in Table |l[ From the four independent transformation pa- 
rameters, only three combinations survive: e + e, 6 + 6, and 6 — 6. Therefore, 
the coefficients of 'ilìD^ip, '{['{■y ■ D)'^'ip, and '4>[yoDQ,'y ■ D]^p can be chosen arbi- 
trarily. The last of these has no redeeming features, so 6 — 6 should be chosen 
so that it never appears. 

The other two redundant interactions are again used to solve the doubling 
problem. The Dq terra is used to eliminate states that would make contributions 
to the fermion propagator proportional to (—1)*; the factors 1 it 70 in eqs. ( |1.6D 
or ( |2.3D provide the unique choice. Low-energy states with pi ~ vr/a are lifted 
by adding the interaction proportional to r<j in eqs. (|l.6| ) or (|2.3| ). When the 
mass is nonzero, it may prove convenient choose to be a function of mga, so 
we leave it arbitrary. 

As with axis-interchange symmetry, the chromomagnetic and chromoelectric 
interactions are not redundant. Their couplings can be used to remove cutoff 
effects once the doubling problem has been eliminated. 

Thus, without axis-interchange symmetry there are eight interactions up 
to dimension five, one of which goes with wavefunction normalization (e.g. 
tpjoDoìp). Three couplings are redundant; two can be used to solve species 
doubling {ìpDQÌp and ìpi^j ■ D^ìp), and the other to eliminate 'ip[yoDQ,j ■ D]ìjj. 
The other couplings are fixed by the fermion mass (nioipip) and three physical 
improvement conditions {Cip^-Dip, icBip'^-B'ip, and CE^ot-Eip)- 

Redundant combinations of higher-dimension interactions can be exposed 
by generalizing the transformation of eq. ( ^.21 ). In particular, after dispensing 
with axis-interchange symmetry, it is possible to transform away interactions 
with higher time derivatives of ip and ■ì/', in favor of spatial derivatives of ip 
and B and E, and time derivatives of the latter. Indeed, any action with 
the Wilson time difference — the first line of eq. — has an easy-to-construct 
transfer matrix. (This is reviewed in sect. |5[) Consequently, it is possible to 
implement eq. ( |1.3| ) by adding "spatial-only" interactions to Sq. 
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4 On-shell Correlation Functions 



We now turn to the mass dependence of the tree-level couphngs, generically 
denoted Cn^(moa) in eq. ([L.4|), needed to bring the action closer to the renormal- 
ized trajectory. This section uses the fermion propagato! to obtain the relation 
between the physical mass and the coupUng moa, the correct tuning of the cou- 
pUng C(moa), and the normahzation of the field tp{x). Since we are interested in 
the full mass dependence, we do not expand in the fermion mass. Sect. ^ uses 
the Hamiltonian of the lattice theory to clarify and extend the analysis to cs 
and CE- 

A well-known procedure for determining the couplings [ Q is to calculate 
n-point correlation functions and expand in momentum p. In gauge theories, 
however, it is not known whether lattice artifacts can be removed systematically 
from Green functions off mass shell. Hence, one expands "on-shell" quantities 



instead [ 19|. The (lattice) mass shell specifies the energy at given spatial mo- 
mentum p, so on-shell improvement amounts to an expansion in pa. Previous 
analyses [ H], ^, also expanded in the coupling moa. We simply avoid the latter 
expansion, and thus obtain the full mass dependence. 

The simplest on-shell correlation function is the fermion propagator as a 
function of time and spatial momentum. It is used to relate the bare mass to 
a physical mass and to derive the mass dependence of (. In the language of 
sect. ^, it probes the interactions 'ij^ip and -{jj-y-Dil;, relative to ^/'7o-DoV'- 

Define C{t,p) through 

{^it',p')i;it,p)) = i27rf6{p'-p)Cit'-t,p), (4.1) 
wheie ip{t,p) = a'^ X^x ^~*''^V'(^; ^) and similarly for ^/5(t, p). Then from eq. ( ^.31) 

TT 27r sin Po + iCj-S -|- mo + 1 - cospo + ^i^sCP^ 

where Si = a^^sinpjO and pi = 2a^^ sin(pja/2), but for brevity eq. ( |4.2| ) is 
given in lattice units. To integrate over pQ, proceed as follows: rationalize the 
denominator; for t > let z = e*^°, and for t < let z = e~^^'', yielding a 
contour integrai over the circle \z\ = 1; apply the residue theorem to obtain 

r<i4. \ -7 -EìtHo sìgntsìnhE -iCfS + mo + 1 - cosh E + ^rsCp^ , . 
^^'^P) = ^^' ^^^[^ (4.3) 

for t 7^ 0,^ where (restoring a) 

„shi.a=l+<!=2^±M^!^2W^ (4.4) 

2(1 + moa + \rsCp^a^) 



Cit,p)=l ^— , _ ^ (4.2) 



To obtain C(0) from eq. (^j), replace 70 signt by 1 on the right-hand side. 
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implicitly defines the energy of a state with momentum p. The residue Z2 (p) is 
given below in eq. ( [4.12 ). 



Expanding the energy-momentum relation in powers of pa yields 

E' = M? + ^y + ..., (4.5) 



where the "rest mass" 
and the "kinetic mass" 



Mi=E{0), (4.6) 

= id^E/dpX=o- (4.7) 

(Any axis i will do to define M2, by spatial axis-interchange symmetry.) The 
relativistic mass sheh has rug = Mi = M2, and it terminates at p^. From the 
tree-level eq. (4.4) 

Mi = a"Mog(l + moa), (4.8) 



and 



2e ^^v^^ (^^gj 



moa{2 + mga) 1 + mga 

Eq. (f4.8D shows how to adjust moa so that = Mi. Similarly, eq. ( |4.9| ) shows 
how to adjust C and so that rug = Setting Mi = M2 and solving for C 

yields the (tree-level) condition (setting o = 1 again) 



// r^mo(2 + mo) y mo(2 + mo) r.,mo(2 + mo) 
^~VV 4(l + mo) y 21og(l + mo) 4(1 + mo) ' ^' ^ 

The dimension-five coupling is treated here as redundant; it is determined 
not by physics, but to solve the doubling problem. To alleviate doubling any 
rs(0) > will do, and the most naturai choice is rs{0) = 1. 
For small mass the Taylor expansion of eq. (|4.10D is 

C = 1 + 1(1 - r,(0))mo - ^[1 - 3r,(0)(2 + r,(0)) + Ur'MWo + Oiml). (4.11) 

At niQ = the redundant coupling drops out, leaving C(O) = 1 unambiguously. 
On the other hand, the full mass dependence of C can only be prescribed hand- 
in-hand with r^. The origin of the link between the two couplings is that both 
the kinetic {'ipyD'il:) and Wilson (^A^^^i/') terms contribute to E"^ at O(p^). This 
and analogous links between couplings' mass dependence are examined further 
in sect. ^ and Appendix 

Beyond tree level (in perturbation theory or in Monte Carlo calculations) 
one would fune C according to the same physical principle that led to eq. ( |4.1C| ): 
determine the momentum dependence of the energy of a suitable state and 
demand that Mi = £^(0) and M2 = {d^E/dpf)'^^ he equal. 

^In the Kt-Ks parametrization (eliminate mo with eq. ( |2.2| )) this condition is an imphcit 
transcendental equation. In the mo-C parametrization one can solve for ( explicitly. 
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When C, and niQa have been adjusted so that Mi = M2 = rug, one can rewrite 



eq. ( |4.5D as E = \Jrn^ + (J-Eiat- Expanding eq. ( ^ ) to p^, one finds the 
lattice artifact 5E\^i ~ p^a^ at smah mass and 5E\^-^ ~ p^a/Ml at large mass. 
To reduce 5Siat further, one must incorporate higher-dimension interactions into 
the analysis. 

Finally, let us return to the residue Z2 in eq. ( [4.31) . In general, the residue 
is a scalar function of four-momentum p, evaluated on shell. With a Euclidean 
invariant cutoff, scalar functions can depend only onp^; on shell, withp'^ = —rn^, 
the spatial momentum p drops out. With the lattice cutoff, however, the mass 
shell is distorted, cf. eq. (|4.4| ), so three- momentum p dependence can remain. 
Indeed, after integrating eq. (|4.2|) over po one finds 



Z2(p) = (l + moa + ir,Cp'a') (4.12) 

Normally one identifies the residue with a (re)normalization of the fermion field. 
Now, however, it is appropriate to expand Z2{p) = Z2 + O(p^), where 

Z2 = (1 + moa)-i = e-*^^^ (4.13) 

Then '^^'^il){x) = e^^°'/'^ilj{x) has the canonical normalization. In the hopping- 
parameter notation the canonically normalized field is ^Jl — QrgUs ipn- This nota- 
tion shows clearly that the approach to the static limit, Ks <^ 1, smooth. Indeed, 



eq. ( 4.13 ) captures the dominant mass dependence of the field normalization to 



ali orders in perturbation theory, cf. sect. S and ref. [ 18|. 



One might ask what to make of the momentum dependence of Z2, when the 
action is improved to higher dimensions. The residue itself is not observable; 
physical quantities are given by ratios of n-point functions and the propagator. 
With the correct on-shell improvement, the p^ dependence of untruncated n- 
point functions combines with that of Z2 to yield the desired results (to the 
order considered). 

5 The Hamiltonian 

This section introduces another method for deriving conditions on the couplings 
in the action. The strategy is to obtain an expansion in the lattice spacing for 
the Hamiltonian. For concreteness, we focus on the action S = Sq + Sb + Se- 
The couplings are then adjusted so that the Hamiltonian of the lattice theory 
is equivalent to the Dirac Hamiltonian. The idea is conceptually the same as 



on-shell improvement, because the "spectral quantities" of ref. [ 19] are just 
eigenvalues of the Hamiltonian. But since the Hamiltonian is an operator, it 
contains the Information of infinitely many quantities, rather than the finite 
number accessible when one computes correlation functions. 

This approach reproduces the Cn{mQa) derived with on-shell correlation func- 
tions. But the analysis is explicitly relativistic, if noncovariant, so one sees 
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clearly that the resulta are general. On the other hand, we have not attempted 
to extend the method to four-fermion operators, or to higher orders in gQ. The 
calculations required by those extensions seem simpler with Feynman diagrams. 

There is a further conceptual advantage to the Hamiltonian. Lattice field 
theories are almost always formulated in imaginary time. The interpretation of 
the results in real time hinges on a good Hamiltonian fixing the dynamics of 
the Hilbert space of states [|^. Hence the implicit, but seldom stated, goal of 
improvement is an improved Hamiltonian; this section merely takes direct aim 
on that goal. Moreover, once one accepts the centrai role of the Hamiltonian, 
one appreciates why a satisfactory Hamiltonian H implies a satisfactory time 
evolution e~^", no matter how large Ha is. 

In lattice field theory the Hamiltonian is defined through the time evolution 



operator, or "transfer matrix" [|2[. Therefore, sect. |5.l| starts by reviewing and 
extending the construction of ref. [ ^ to the actions Sq and S = Sq + Sb + Se- 
A by-product of this analysis is the demonstration that there is no need to 
improve the temporal derivative in eq. ( |2.3D . This feature is familiar from the 
static and nonrelativistic formulations. It is a special blessing here, because a 
temporal next-nearest-neighbor interaction would introduce unphysical states [ 
0], and at large m^a the physical and unphysical levels cross. With the transfer 



matrix in hand, sects. |5.2| and |5.3| develop an expansion in a for the Hamiltonian 
itself. 

5.1 Construction of the transfer matrix 

The transfer-matrix construction with two hopping parameters differs little from 
the usuai case [ The transfer matrix acts as an integrai operator in the 
space of gauge fields; in the C/q = 1 axial gauge a wave functional ^t{U) at time 
t evolves to 

^t+l{V) = j \{dUn,^ÌC{V,UWU) (5.1) 

n,i 

at time t + 1. The wave functional 0,t{U) is also a vector in the fermion Hilbert 
space. For the standard gauge action the kernel may be written 

ICiV, U) = fXv)TMv)ICG{V, U)TG{U)fF{U). (5.2) 

The factors arising from the fermion action Tp^ are operators in the fermion 
Hilbert space. The factors arising from the gauge field, Tq and /Cg, are given in 
ref. [ 24 1; in the following, they do not play a cruciai role, so we do not discuss 
them further .0 



^Different from ref. [ ^ is the convention for the factors (1 ± 70) in the action (compare 
eq. (1.6) with eq. (2) of ref. [ ^). With our convention it is naturai for time-ordering to place 
later times to the left. Thus, the kernel KiV, U) transfers the field from U at time t to V aX 
time t 
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The fermion operator for action Sq can be written 



f^(U) = e-^Hf^)e-5^o(^) det{2KtBu)'^^ (5.3) 
where (cf. ref. [|2|]) 

Ho{U) = ^Mu^, (5.4) 

Hj{U) = C^^{l-^oh-Du^, (5.5) 

in a matrix notation in which ^ and ^ are vectors and Bfj, Du, and À4u are 
matrices depending on gauge field U. The vectors and matrices of this notation 
are labeled by spin, color, and space. The covariant difference operator D is as 
in eq. ( |2.5D and 

B=l-rsKsY,iTi+T^Ù, (5.6) 



e 



M 



B 



l + mo-ir,CA(3). (5.7) 
The operators ^ and ^ = ^^70 obey canonical anti-commutation relations 

{^™a,^L} = Wafe, (5.8) 

where m and n label spatial sites and a and ò are multi-indices for spin and 
color. The fields corresponding to these operators are related to the originai 
fields by 

^ma = nbV'nfe. (5.9) 

This discrepancy in normalization between the integration variables in the func- 
tional integrai and the canonical operators in Hilbert space demonstrates again 
that the normalization convention for the field ^(x), cf. eq. ( |2.1[) , is arbitrary. 
On the other hand, the propagator of ^' has unii residue at tree level, and a 
perturbative series ^2^1 = 1 + gQZ]^^{mQa) + • • • beyond tree level. 



The generalization of eqs. (|5.2| )-(5.9) to include the chromoelectric interac- 



tions suffers from a technical difficulty. Usually one uses the "four-leaf clovers" 
in eq. ( ^.12 ) as the lattice approximants to the chromomagnetic and chromoelec- 



tric fields. For the chromomagnetic interaction, this choice poses no problem, 
because B involves link variables from one timeslice only. For the chromoelec- 
tric interaction, however, the time-space four-leaf dover involves link variables 
from three timeslices. In that case, the construction of the gauge-field transfer 
matrix is more complicated, and, if the improved gauge action is any indication, 
it may no longer be positive [ |2^ . 

To avoid this complication one can define a chromoelectric field on only two 
timeslices. Consider 



Se2 = -CEl^s^'4'n,t + ^o)o:-^^n,t-l/2 + è(l ~ ^0)"-^^™,t+l/2 Ì^n,t, 

(5.10) 



n,t 
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where ^ 

En,t±l/2;i = ±^ XI SÌgn(2:) T^qTiT^qT^^ - h.c, (5.11) 



is defined on a tooleaf dover. The projection operators ^(1 it 70) in eq. ( 5.10 ) 
are chosen by analogy with the Wilson time derivative, cf. eq. ( |2.3D , and as a 
result the standard transfer-matrix construction goes through with minor mod- 
ifications. The two-leaf version Se2 differs from the four-leaf version Se by an 
interaction of dimension six, so it should not alter the tree-level tuning ol ce- 

Extending the transfer-matrix construction to 5o+S's+5£;2 (eqs. ( |1.6| ), ( |1.7D , 
and (|5l0| )), one finds the following changes. The chromomagnetic interaction 
modifies the matrices B and M to 

B = l- rsKs X(?i + r_i) - icBKs^-B, (5.12) 

i 

= ^ = 1 + mo - K (r-sA(3) + ìcb^S-b) . (5.13) 

Except for the new B, eqs. ( |5.8D and ( |5.9D stili hold. The chromoelectric inter- 
action, eq. ( ^.10 ), modifies the fermion operator Tp so that it depends on initial 



and final gauge fields U and V: 

TeÌV, U) = e-^^(^.^)e-5^o(^) àei{2KtBufl' (5.14) 
with Hq as in eq. ( p.4[ ) and M. from eq. ( |5.13| ), but 

Hi{V, U) = Cl'i(l - l^){l-Du - \cECx-Ev,u)^, 
Hj{V, U) = CÌ>\{1 + 7o)(7--Dy - \cEOc-Ev,u)Ì!, 



(5.15) 



where the subscripts on D and E specify the spallai link fields, out of which 
they are constructed. The sign of the chromoelectric term in h\ can be checked 
as follows: in our sign and i conventions t'^'^ = —f^ and Euclidean electric fields 

are anti-Hermitian operators in the gauge- field Hilbert space, È""^ = —È"" 



Comparing eqs. ( |5.13|) and (|5.15| ) with eqs. (5^) and (|5.5|) , respectively, one 



notices a pattern emerging. Interactions with block-diagonal Dirac matrices ap- 
pend to e-'^, whereas those with block-off-diagonal Dirac matrices modify Hi. 
This pattern depends only on the special Wilson time derivative and the tech- 
nical assumption that ali interactions live on only one or two timeslices. It 



proves that there is no need to alter the temporal derivative in eq. (2.3): higher- 
dimension "spatial" interactions are enough to achieve on-shell improvement, as 
asserted at the end of sect. 0. 
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5.2 Small a expansion (general considerations) 

From the transfer matrix one would like to derive the exact lattice Hamilto- 
nian H = — log/C. Of course, with this defìnition the Hamiltonian cannot be 
represented by a finite number of locai operators. According to the Symanzik 
philosophy, however, one ought to expand it in powers of the lattice spacing. 
After obtaining the transfer matrix, (higher) time derivatives are no longer a 
concern, so the lattice-spacing expansion will hold if the quantities Da, Ba^, 
Ea?, . . . , are small. 

One can anticipate the expansion by enumerating the terms allowed by sym- 
metry: 



)2 

cont 



H = ^ bo{moa)mq + 61 (77100)7-1? cont + ab2{moa)D^ 

+iabB{moa)'S-B cont + abEÌmoa)a- E cont H ^, 



(5.16) 



where the subscript "cont" refers to an underlying continuum gauge field; below 
we usually suppress this subscript, for brevity. The coefìicients òj depend on moa, 



and since eq. ( 5.16 ) is to be interpreted as an expansion in a (rather than 1/mo), 



the òj for small mQU must be O((moa)^), with p nonnegative. The coefficients 



for the action Sq + Sb + Se, given in sect. 5.3, satisfy this requirement. 

The general objective is to adjust the couplings so that eq. ( 5.16 ) takes the 
relativistic Dirac form, i.e. Òq = òi = 1 and b2 = bs = bE = ■ ■ ■ = 0- But, based 
on the considerations of sect. ^, there must be some leeway in the redundant 
directions. In the operator formalism adopted here, unitary changes of variables 
are possible, and these play the role of the isospectral transformation, eq. ( |3.2D . 



Under a change of variables the Hamiltonian becomes 

È' =Ù{Ù + dt)Ù-^, (5.17) 

where dt is a derivative with respect to imaginary time, and lA is the unitary 
operator implementing the change of variables in Hilbert space. Consider, for 
example, the following transformation: 

^ ^ exp(-<i7-i:))^', 

(5.18) 

^ ^ ^'exp(-a^i7-I?), 

for which 

= exp (a^i^^-D^') . (5.19) 

Such transformations are familiar from studies of the nonrelativistic limit of the 
Dirac equation, where they are called Foldy-Wouthuysen-Tani transformations [ 
pO| . Their characteristic feature is that the exponent is always a block-off- 
diagonal Dirac matrix. 
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The transformed Hamiltonian H' has an expansion of the same form as in 



eq. ( 5.16 ), but with transformed coefficients: 

= ^0, 

ò'i = hi- 2mga6oCi, 

ò'2 = 62 - 2616 + 2mqabQÌl, (5-20) 
b'B = bB- 2616 + 2mgaòoCi, 
b'E = bE-Ù- 

In light of the transformations, it is, therefore, enough to adjust moa, (, r^, cb, 
and Ce, so that for some (hidden) vaine of ^1 the transformed Hamiltonian takes 
the Dirac form H' = ^{niq + ■y ■ D)'^ . That means that one wants 6q = b'^ = 1 
and b'2 = b'B = b'j^ = ■■■ = 0. 

The Foldy-Wouthuysen-Tani parameter ^1 drops out of on-shell quantities. 
It is preferable, therefore, to parameterize the redundant direction by one of the 
couplings. To this end, it is efìicient to note that the following combinations of 
the òj's do not depend on ^1: 

Bo = bo = ò'o, 

Bi=bl- 2mgabob2 = b'^ - IrUqab'^b'^, (5.21) 
BB = b2-bB = b'2 - b'^. 

A Hamiltonian unitarily equivalent to the Dirac Hamiltonian is then obtained 
whenever 



Eqs. (5.21) do not contain an invariant corresponding to bs- This is analo- 
gous to the result, eq. ( ^.11| ), that the general mass dependence of C, can only 
be determined hand-in-hand with Vg- In the present language, that connec- 
tion arises as follows. Consider truncating eq. ( ^.16| ) at dimension four. Then 
only òo and òi remain. The Foldy-Wouthuysen-Tani transformation is 0(pa), 
and superficially not worth considering. If one introduces it anyway, one sees 



immediately that 61 changes, and in the transformation law, eq. ( 5.20 ), one 
power of a has combined with the fermion mass to give niqa. At m^a = the 
condition 61 = 1 is enough to determine ("(0) unambiguously. But to obtain fully 
the mass dependence C('^o) one must consider simultaneously the Inter actions 
^^■D^ and ^{-y-D)"^^. 

A similar fate awaits the coefhcient ò_b and its coupling ce- Consider a 
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two-parameter Foldy-Wouthuysen-Tani transformation 



^ ^ exp{-a(,i'-f-D - a^^EOL-E)-^, 

^ ^ ^'exp(-a^i7-i:) - a?£,EOL-E). 

The new parameter introduces changes that are superficially O(p^a^) 
other coefficients are unaffected by but 



(5.23) 



The 



ò'e = hE-Ìi- 2m„o6oCE- 



(5.24) 



Again, one power of a has combined with the fermion mass to give m^a. Thus, 
the condition hE — Ìi = is enough to determine only C£;(0). The full mass 
dependence of ce can only be revealed by considering simultaneously ^a-E"^ 
and the dimension-six interaction q:-£^}^. This analysis is deferred to 

Appendix 

The next subsection adjusts moa, fs-, and cb to ensure eq. ( ^.22 ), and ce{^) 
to ensure 6'^(0) = 0. For heavy-light systems, the resulting lattice theory has 
cutoff artifacts of 0(AQ(-;£,a^) and, only when mqa ^ 1, (yet smaller) artifacts 
of 0(Aq(jj-) a/mg) and of 0(AQQj-,/mQ) as well. See sect. |6|for details. Moreover, 
for quarkonia, the lattice theory is similarly correct through 0{v^). 

5.3 Small a expansion for Sq + Sb + Se 



Combining eqs. ( ^.2| ) and ( 5.14| ), and omitting factors that depend only on the 
gauge field, the fermion Hamiltonian of the lattice theory is 



H 



log 



^-\Ho{V)^-H]{V,U) 



(5.25) 



where Hq is specified by eqs. ( |5.4| ) and ( 5.13| ), and Hi is specified by eq. ( 5.15| ). 
To derive an expression for the fermion Hamiltonian, one must coalesce the four 



exponents in eq. (5.25) into one. Owing to nontrivial commutators between Hq, 
Hi, and h\, this is too difhcult in general. But through order D'^a? , Appendix |^ 
achieves the desired result by a trick. There the field theory is mimicked by a 
toy model with the same algebraic structure but only two degrees of freedom. 
In the toy model one needs only to take the logarithm of a two-by-two transfer 
matrix, and expand the result in powers of a. 
For small Da the Hamiltonian becomes 



H 



Mi 



C 



(rsA(3) +ìcbT,-B 



2(1 + mo) 

Kh{mo)Q - CV2(mo)e2l Ì! + 0(p3a2), 



where 



e = i(7-£>cont + è(l-C£;)a--E 



cont j 



(5.26) 



(5.27) 
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The rest mass Mi and the terms in parentheses come from expanding M, and 
the functions fi are extracted from the toy model: 



/i(^) 



2(1 + x)log(l + x) 
x(2 + x) 



f2{x) 



1 



21og(l + x) x{2 + x) 



(5.28) 



Note that /i(0) = 1 and /2(0) = i. 

In the spirit of an underlying continuum gauge field one can identify A*-^) 
with i^cont , DxD with -Bcont 1 and [dt , D] with .Econt • With these identifications 
one can cast eq. (5.26) into the form of eq. ( 5.16| ). Thus, the Hamihonian of the 
action Sq + Sb + Se has coefìicients 



bo 
hi 

b2 



Mi /mg, 
C/l("^o), 

ef2{mo) 

C'/2(mo) 



and the invariants Bi are 



2(1 + mo)' 

cbC 
2(1 + mo)' 

= i(l - CE)Cfi{mo), 

Bo = Mi /mg, 

Bi = M1/M2, 

1 1 



(5.29) 



(5.30) 



Bb 



2Mb 2M2 
The masses Mi and M2 are as before, and 

1 , cbC 



+ 



(5.31) 



Mb mo(2 + mo) ' 1 + mo 

After imposing eq. ( |5.22 ), Mi, M2, and Mb ah equal the physical mass. 

The mass dependence of the couplings fohows immediately from eqs. (5.30) 
and ( 5.22 ). The requirement Bi = 1 imphes 



c 



rsmo{2 + mo) 



+ 



mo(2 + mo) rsmo(2 + mo) 



4(1 + mo) / 21og(l + mo) 4(1 + mo) 
precisely as in eq. (|4.1C| ). The requirement Bb = imphes 

CB = Ts- 

Finahy, for small mass the chromoelectric couphng should be tuned to 

CE(0) = i(l+r,(0)), 



(5.32) 

(5.33) 
(5.34) 
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to enforce b'^{0) = 0. With the axis-interchange invariant boundary condition 
rs(0) = 1, one thus recovers the action of ref. [|6|, with ({0) = rs(0) = ce{0) = 
cb{0) = 1. 

Our analysis has not yet specified the relevant couphngs gQ and m-oa. They, 
of course, are fixed not by theory but by experiment. In the Hamiltonian lan- 
guage, the bare mass moa is adjusted so that Bq = 1, i.e. Mi = niq. Then 
the improvement conditions, eqs. ( ^■32D -( ^.33D , guarantee that iriq = M2 = Mb 
also. 

There is a special case of eqs. ( 5.32| )-( 5.34 ) that is of at least passing interest, 



namely the one for which the Foldy-Wouthuysen-Tani parameter = 0. This 
is obtained by choosing so that (untransformed) 62 = 0: 



2(1 + mo) 



mo(2 + mo) log(l + mo)' 
Then the condition òi = 1 requires 



(5.35) 



r = ?7ìo(2 + mo) 

2(l + mo)log(l + mo)' ^ ' ^ 

and the condition 6^; = requires 

CE = 1, (5.37) 

and, as before, the condition 6^ = requires cb = rg. After substituting 
eq. ( |5.35|) into eq. ( |5.32| ) one re-obtains the right-hand side of eq. ( ^.36 ). Ap- 



pendix ^ shows that — with and ( from eqs. ( 5.35| )~( [5.36| ) — ce = 1 can be 
maintained for arbitrary moa. 



6 Truncation Criteria Revisited 

This section reexamines criteria for truncating a cutoff theory, with some empha- 
sis on the errors left over after truncation. The analysis of the previous sections 
takes the scaling dimension of the interaction as a guide. For massless quarks 
that is certainly correct. But the most appropriate organization may vary when 
the same cutoff theory is applied to different physical systems. Thus, conclusions 
about the accuracy of a massive- fermion action must be refined, after deciding 
whether the action is to be applied to heavy-light systems or to quarkonia. 

After the couplings have been adjusted to some practical accuracy, the 
Hamiltonian (possibly after a Foldy-Wouthuysen-Tani transformation) is 

iJ = l'(mg + 7oAo+7-£>)^ + 5Aat; (6.1) 

the Coulomb potential appears if one transforms to a gauge without Aq = 0. 
A lattice artifact dHi^t remains, because one cannot exactly incorporate infinitely 
many terms into eq. (|1.3|). 
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One can estimate the errors induced by by treating it as a perturbation. 
There is an advantage to estimating cutoff effects from the Hamiltonian. In the 
action formahsm, eq. ( |1.5| ), it may not be clear how the time discretization 
trickles down to physical quantities. But by proceeding through the transfer 
matrix these effects are treated exactly. 

From the hne of argument leading to eq. ( 5.16| ), one expects that J-ffiat 



consists of operators multiphed by mass-dependent coefficients 

oo 

(^i?iat = E«'" E 9oh^Kmoa)Hr., (6.2) 

ri l=Ln+l 

where the power Sn = dim Hn — 4, and L„ is the number of loops ah'eady under 
control. One can determine the effect of SHia,t on a physical quantity from 
order-of-magnitude estimates for the operators Hn and general properties of the 
coefficients (moa). While the former depend on the physical process under 
study, the latter are process independent. 

The dimension-five, tree-level coefficients have two important properties, 
which we believe are generic. First, at asymptotically large moa the tree-level 
coefficients either approach a Constant or fall as a power of l/(?7ioa). An analysis 
of higher-order Feynman diagrams (sect. ^) shows that tree-level patterns persist 
to ali orders in perturbation theory. Indeed, the asymptotic behavior is presum- 
ably a consequence of the heavy-quark symmetries obeyed by ali lattice actions 
under consideration. Second, the coefficients always contain the recurring ingre- 
dients 1 + mga, moa(2+moa), and log(f +moa) in a way that makes implausible 
any combination that would blow up at an intermediate value of moa. Indeed, 
ali evidence suggests that the functions h{mQa) are smaller than their low-order 
Taylor expansions, once moa^ 1. 

Let US now discuss the typical size of the operators in the Hamiltonian. 
Table | gives ballpark estimates for the dimension-three, -four, and -five inter- 
actions for three systems: those in which ali quarks are light, those with one 
heavy quark, and quarkonia. The row labeled Eq in Table || gives the non- 
trivial dynamical scales, to which artifacts should be compared. In all-light and 
heavy-light systems, the estimates start from naive dimensionai analysis, but 
heavy-quark bilinears with an off-diagonal Dirac matrices are Aqcd/^tiq times 



smaller stili. In quarkonia, the estimates are those of ref. [ 12|, with v denoting 
the typical velocity of the heavy (anti-)quark in the bound state {v ~ as^mq)). 

A conservative estimate of the artifact is then as follows: Choose a system, 
multiply by a'^", and compare to £"0. The coefhcient ò(moa) is a number of 
order 1 (or less), for any value of moa, so its numerical value does not affect 
the (conservative) conclusion. If, after suitable adjustment of the couplings, one 
finds b ~ TTigO for m^o <C 1, or ò ~ l/(mQa) for mga ^ 1, the artifact might be 
even smaller. 

Consider the chromoelectric interaction as an example. For the sake of argu- 
ment, suppose that the rest mass Mi and the kinetic mass M2, and hence mga 
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Table 2: Estimates of the size of dimension-three, -four, and -five interactions 
in systems with only light quarks, with one heavy quark, and in quarkonia. The 
latter two columns use mg to emphasize the heavy-quark mass. 



Hn 


only hght 


heavy-hght 


quarkonia 




Aqcd 


Aqcd 






1 


1 


1 




A-QCD 

^^QCD 

^^QCD 
A2 

^^QCD 


^^QCD 
A2 

^^QCD 


2 2 
TTIqV 

2 4 

2 4 
TUqV 



and have been adjusted nonperturbatively. li ce is not adjusted correctly, 
then the (transformed) coefficient h'^ of the chromoelectric terni in the Hamilto- 
nian does not vanish. Then, relative to the corresponding Eq, there are artifacts 
of ©(Aqcdo) for all-light, O^K^Q-Qa/mq) for heavy-light, and ©(mgaw^) for 
heavy-heavy. If instead C£;(moa) is adjusted to ce{^) in eq. ( |5.34 ), the artifacts 
in all-light systems fall to O(mgAQCDfl^)- With heavy quarks the estimates 
depend on mqa. If a is so tiny that mqa <C 1, then b'^ ~ n^-QO,, and the chromo- 
electric artifact is reduced to 0(Aq(2;£,o^) for heavy-light and to O(mQa^v^) for 
heavy-heavy. But if mga ^ 1, it turns out that b'^ either remains Constant or falls 
as l/(mQa), depending on the mass dependence of the redundant coupling r^. 
The artifacts are then either ©(Aq^pa/mQ) for heavy-light and O(mQav^) for 
heavy-heavy, or l/^rriga) times smaller. 

The appearance of l/(mQa) in coefìicients, in addition to the Aqcd/^q in 
heavy-light dynamics, makes the error analysis of heavy-light systems somewhat 
delicate. Since the l/{mQa) behavior arises only if mga ^ 1, it leads only to 
errors that are smaller than the usuai discretization errors, relative to Eq. On 
the other hand, occasionally one is interested in effects that are subleading in 
the heavy-quark expansion. For a given lattice action, such quantities may have 
a larger relative error. For example, even with ce{0) adjusted correctly, the fine 
structure of the heavy-light spectrum, which is 0(Aqqj-, /mg), suffers a relative 

error of order (A^cDa/"^Q)/(^QCD/"^Q) ~ ^QC^a. 

Similar comments apply to quarkonia. Though the chromomagnetic and 
chromoelectric interactions are of order itiqv'^ /niQV^ ~ v"^ smaller than the spin- 
independent kinetic energy, they introduce relative errors on spin-dependent 
structure of order mQV^/mQv'^ ~ 1. A full 0(^;'^) analysis requires a few 
dimension-six and -seven interactions, which we consider in Appendix 

Once the dimension-five couplings cb and ce have been properly adjusted, 
lattice artifacts remain from dimension six and higher. Table |3| lists bilinear 
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Table 3: Estimates of the size of dimension-six bilinear interactions in systems 
with only light quarks, with one heavy quark, and in quarkonia. For the actions 
considered here, the interactions below the hne do not arise at tree level. 



Hn 


only hght 


heavy-hght 


quarkonia 




^^QCD 


A3 

^^QCD 


S 4 
ITIqV 




^^QCD 








A3 

^^QCD 


^qcdMq 


S 4 
ITIqV 




A3 

^^QCD 


^qcdMq 


mqV^ 


i^'7oS;-(L>oS)* 


A3 

^^QCD 


^qcdMq 




^'7-(L'o-E)^' 


^QCD 


^qcd/''^q 


niQV^ 


^-fo{D-E - E-D)"^ 
i^-y-{D X B + D X B)"^ 


^QCD 
^QCD 


^QCD 





Table 4: Estimates of the size of dimension-six four-fermion interactions in sys- 
tems with only light quarks, with one heavy quark, and in quarkonia. The role 
of these interactions can only be treated consistently in concert with the gauge- 
field action. The interactions above the line can then be considered redundant, 
while those below the line do not arise at tree level. 



Hn 


only light 


heavy-light quarkonia 




A3 

^^QCD 


^QCD rn^v^ 




A3 

^^QCD 


^QCD ^qV^ 


{^l;^t-^/J)^ E / 7/. 


^QCD 


^QCD ^qV^ 



operators that can appear in the Hamiltonian. The conservative estimate of 
the absolute errors caused by these operators is to multiply Table ^ by a^. 
When niQU ^ 1, however, some of the contributions may be, as before, a factor 
of l/{mQa) or l/(mQa)^ smaller. But, again, this subtlety is only cruciai when 
quantities subleading in the heavy-quark expansion are at issue. 

The four-fermion interactions, listed in Table ^ are also of dimension six. 
To generalize the analysis of sect. ^ to encompass these operators, one must 
simultaneously treat dimension-six gauge-field interactions [ ^. The result is 
that {'tp^ot'^tp)'^ and (■07jt°'0)^ are redundant [|^, even without axis-interchange 
symmetry. The other four-fermion interactions arise first at the one-loop level. 

Let US summarize the main points of this section for heavy-light spectroscopy 
with action Sq + Sb + Se- After the tree-level adjustments of sect. 5^ have been 
applied, the largest remaining lattice artifacts are ©(^gAgcDo) from the one-loop 
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maladjustment of Sb + Se and 0(AQQj-,a^) from unadjusted dimension-six inter- 
actions. The mass dependence of the artifacts is solely in the coefìicients h{mqa), 
which is a number of order unity at any mass. 



7 Electroweak Perturbations 

This section extends the formahsm of the previous sections to the two- and four- 
quark operators of the electroweak Hamiltonian, which may be treated as a first- 
order perturbation to QCD. The construction of the renormahzed (or continuum- 
hmit) operator is analogous to the construction of the renormahzed trajectory. 
Let O denote the continuum operator. Then 



O = Zo{{moa},gl) 



Oo + X! Cn{{mQa},gl) 0„ 



(7.1) 



where the sum runs over ali lattice operators 0„ with the same quantum numbers 
as C[[] Like the couplings in the action, the coefficients Zq and C„ are functions 
of the relevant couplings, ali fermion masses {mg a} and the gauge coupling Qq. 



Eq. (7T) is general, but we again consider perturbative expansions in Qq, 

oo 

Zo{{mQa},gl) = ^5g'zgl({moa}), 



(=0 

oo 

Cn{{mQa},gl) = ^gl^cl\{moa}) 



(7.2) 



and focus on tree level. Previous work [0 applied to small masses, but we treat 
the mass dependence of Z^'({moa}) and c'Ì}{{moa}) exactly. We also do not 
impose axis-interchange invariance in classifying the lattice operators On- 

The coefìicients Zq and Cn can be determined from low-momentum matrix 
elements of ali 0„, analogously to sect. |^. In perturbation theory it is enough 
to compute matrix elements between quark, anti-quark, and gluon states, both 
with lattice and MS regulators. It is essential to impose consistent normalization 
conditions. Appendix ^ derives external-state rules for lattice perturbation 
theory. There one finds that the contraction of ip{x) with a normalized fermion 
state corresponds to a factor u\a,t{^,p)M{p)., where «lat is a normalized spinor 
on the lattice mass shell. The factor 

//u(p) -cosh^;^ 



where (for Sq) fj,{p) = 1 + niQa + \rsC,p^a?. A relativistic theory has instead 



Ureìi^,p)ynT'q/E, where «rei and E comply with the relativistic mass shell. 

^By convention, the zeroth lattice operator Oo has the same dimension as the continuum- 
limit operator O. The role of the other C„0„ is to remove terms suppressed (or enhanced!) by 
a power of a. The role of Zo is to convert to a preferred renormalization convention. 
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Consider the bilinear operator J-^^ that creates flavor / and annihilates 
flavor g with spin coupling V. At tree level its matrix elements should be 



{q' (e, P') I kne, P) ) = <i {i' , P')^Kei (e, P) ^fn^jn^ jE^, ó'^a , 
{r{e,p')\j/'\qHC,P)) = -vU^,p)rv'^,i{e,p')Vmamb/E'^Eb 6'f6-9^ 
m/'\q''{tp)^{e,p'))=vU^\p')ruUC,p)^mam,/EaE',ó'f6'^^, 
{q^iC, p)q\Ì' ,p')\JI'\^) = p)Tvl,{i',p')Jmamb/EaE', S^f 5^9 , 



(7.4) 



where e'^^ is the energy of flavor / with momentum p^'^ and mass mj. Note 
that the relativistic spinors «rei and Urei appear on the right-hand side. 

With the right-hand side of eq (7.4) as a target, we now consider lattice 
operators 0„. The simplest lattice bilinear with the correct dimension and 
quantum numbers is 

Jl3{x) = Ìjf{x)T^3^x), (7.5) 

which corresponds to Oq in eq. ( [Z.!]) . Recali that ìI){x) is the field appearing in 
the mass form of the action, eq. ( |2.3D . At tree level the matrix elements are 



('z'(c',p')i^/'k'^(e,p)) =<t(e',p')r<t(e,p)AAa(p)AAfe(p') s^fs-^^ 

{<fÌe,p')\4'\q'{C,p)) = -vUC,p)rvUC',P'Wa{p'WbÌp) 6'^ 5-9, 
{0\4'\q-{C,p)q'{e,pl)=vUe,p')ruf,,{^,pWa{pWb{p') ó'n-a^ 
{q''itp)q\e,p')\4'\0)=uUC,P)rvUC',P'WaÌp)MbÌp') ó-fS's, 



(7.6) 



where J\ff{p) is the normalization factor of flavor /, cf. eq. ([7.3|). Note that the 
lattice spinors liiat and viat appear on the right-hand side. 

Setting p = p' = 0, the matrix elements differ only because of the fac- 
tors M{0). Thus Z^Jr has the same zero-momentum matrix elements as the 
target J7r, in ali four channels, if the (re)normalization factor 



Zr{mofa,moga) = ^ fif{0)fj,g{0) = exp (i(Mi/a + Miga)^ (7.7) 

This is a tree- level result, but the mass dependence shown here remains dominant 
to ali orders, cf. sect. ^ 

Further terms in the three-momentum expansion cannot be matched with- 
out considering higher-dimension terms in eq. ([Z.!]). At tree level one sees the 



differences between eqs. ( |7.4D and (0) in the factors Àf 7^ ^/mjE and spinors 
itiat ^ lirei- Eq. ( [7.1D can therefore be extended to higher dimension by intro- 
ducing an improved field. To first order in pa consider 

^j(x) = e^^^''l'^\\ ^ adxl-D\i\){x), (7.8) 
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with flavor labels implied. Then 



Jp = ^i(T^f{x) (7.9) 

is the target operator of interest, through first order in pa, if di is adjusted prop- 
erly. Comparing the bracket in eq. with those in eqs. (|a25|) and ( |a27|) , 

one finds 

di = + - (7.10) 

moa(2 + moa) 2M2a 

identifying rug = 

For small mass one finds di oc moa; the only 0(a) improvement needed is the 
normalization factor e^'^'^"'^'^. Ai farge mass, however, the rotation of eq. ( \l.8\i 
becomes important. Analogousfy to the Hamiftonian coefficients discussed in 
sect. ^ when moa » 1, one has di ss l/(2mg). Consequently, the contribution 
of di^-D is essential for computing the l/nig correction to the static limit of 
matrix elements of Jr- Similarly, higher-dimension generalizations of eq. ( [T.Sp 
are needed to obtain l/rrig and corrections of higher order in l/rrig.. 

The improved field ^'i(x) in eq. ( |7.8| ) coincides, through 0(p), with the one 
denoted by ^ in sect. ^. Combining eqs. (|5.9D and ( 5.1S| ), the Foldy-Wouthuysen- 
Tani transformed field is 

^{x) = exp{a^i-fD)e^/^xp{x), (7.11) 

where ^i parameterizes the solution of the tuning conditions. This expression is 
(numericahy) cumbersome, but one may expand consistently the exponentials 
in a. This exercise identifies ^i with di. Indeed, solving b'i = 1 for ^i yields the 
right-hand side of eq. ( [7.10| ), after replacing the rest mass Mi with the kinetic 
mass M2.^ 

The special role of \I' should not be too surprising, because it possesses two 
important properties. First, it satisfies canonical anti-commutation relations 
and is thus properly normalized. Second, its dynamics are given by the Dirac 
Hamiltonian H = ^ (rug + ^qAq + ^■D)'if — at least at tree level and up to O(p^). 
Therefore, any operator built out of the transformed field yields the desired 
matrix elements, also at tree level and up to O(p^). 

Let US conclude this section with some comments on two other Ansàtze for 
the currents. A formai argument based on the Ward identity suggests that a 
conserved current^ is especially suited to the determination of form factors of 
a vector current or the decay amplitude of a vector meson. But although the 
Ward identity implies a certain universality in radiative corrections, it does not 
imply any special mass dependence at tree (or any other) level. 

*The substitution of the kinetic mass M2 for the rest mass Afi is done so that the expression 
remains valid under a nonrelativistic interpretation explained in sect. ^. 

^Both "Noether" and "gauge" currents are conserved; they differ by a^i, terms. 
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With standard Feynman rules and Appendix straightforward algebra 
yields the tree-level on-shell matrix elements. To 0(p) the (conserved) gauge 
current V,9 has matrix elements 



(7.12) 



{o\v^m,p)q{e,p')) = vie, ohM^, o 



sinhMia 



(7.13) 



Mca ' 

through 0{pa), where Mi, M2, and Mb are the tree-level massesj^ but 



1 _ 2C ^ CEC[l + il + mof] ^ ^^^^^^ 



Mg mo(2 + mo) 2(l + mo)2 

The Ward identity asserts that these tree-level masses ali renormalize in a co- 
herent way. But although the "forward-scattering" matrix elements in eq. ( 7.12 ) 
are correct (assuming M2 = Mb = rug), the "annihilation" matrix elements in 
eq. ( 7.13 ) are not (unless m^a <^ 1). We conclude, therefore, that V"^ is not 
useful for determining the decay Constant of a massive vector meson. 
Ref. [ 1^ suggests using a "(four-dimensional) rotated current" 

^r'rot = + è«^)r(l - èa|))V^(x). (7.15) 

To ascertain if j/ matches the target continuum operator, one must evaluate 
matrix elements, as above. The timelike translations in Dq greatly change the 
mass dependence. One finds that Zr,rot<^r rat correctly normalized matrix 
elements only if 

7 (m n rr, ^^ - 4Zr (mp/g, mogo) 

Zr,rot[mofa,moga) - — — . -— — . — r, (T.lfaj 

(2 + smhMi/a)(2 + smhMiga) 

where Zr is the normalization factor of the unrotated bilinear, eq. ([7.7]). More- 
over, when rUga ^ 1 the rotation of eq. ( [7.15| ) must be supplemented à la 
eq. ( |7.8|), w ith the same di as in eq. (7.10). Thus, mass-dependent improvement 
of eq. (|7.15| ) is analogous to improvement of eq. (|7.5|) , but the latter is simpler. 

In summary, the mass dependence of electroweak operators is tractable, if 
one proceeds as follows. First, start with a simple operator Oq and expand its 
on-shell matrix elements in external, spatial momenta small in lattice units. As 
usuai, there is no need to expand in moa. Second, add additional terms CnOn 

^°For the No ethe r curr ent th e terms proportional to cb (imphcitly in 1/Mb) and ce (in 
1/Mg) in eqs. (7.12) and (7.13) would not appear. 
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to correct the momentum dependence of the matrix elements. At least to tree 
level this step can be accomphshed by field rotations, as in eqs. Q and ( |Al7|) . 
Finally, normahze Oq + CnOn to obtain the fuhy renormahzed operator O 
in the desired renormahzation scheme. For example, through 0{pa) the renor- 
mahzed bihnear J'^^ is given by eqs. ([f!s|), (^), and ( |7.10| ). 



8 Beyond Tree Level 

In the previous sections, the niQa dependence of the couphngs in the action is 
derived at tree level. This section considers what happens beyond tree level. 

In perturbation theory the expressions for the masses introduced previously 
become power series in g^. For example, 

oo 

^(0) ^ Mi = m;°i + J2 al'Mf (8.1) 
1=1 



and 



' o2 ITI \ ^ oo 

^""^ ^M, = MM+V,g'MW, (8.2) 



dpf 



p=0 1=1 



where M| and Mg are given by eqs. ([4.8|) and (|4.9| ), respectively. After cal- 
culating the self-energy to l loops, one can extract the coefìicients m|'' as func- 
tions of niQa. The requirement Mi = M2 = subsequently yields the per- 
turbative power series for the couplings and (Based on the arguments of 
sects. m and ^, the Wilson term's coupling should be redundant to ali orders 
in Qq.) In the same vein, the on-shell fermion-gluon vertex function to / loops 

yields Cg(moa) and c^(moa), and electroweak matrix elements to l loops yield 

^^(moa) and él\mQa). 

A complete derivation of one-loop corrections is beyond the scope of this pa- 
per. It is easy, however, to assess two qualitative features: the mass dependence 



of loop diagrams (sect. |8.1|) and the expected size of corrections from tadpole 
diagrams (sect. |8.2D . 



8.1 Mass dependence of loop diagrams 

This subsection shows that the mass dependence of loop diagrams is benign. 
Although we focus on the specific action S = So + Sb + Se, the conclusions hold 
for any action with the Wilson time derivative and arbitrary spatial interactions. 
Actions with next-nearest-neighbor interactions in time are problematic starting 
at tree level [ ^, so they are not considered bere. 

Let US first consider vacuum polarization. At one loop it is easy to see 
that the lattice-regulated Feynman integrals for vacuum polarization are smooth 
functions of the fermion mass. Moreover, for large fermion mass the integrals 
vanish as (moa)~^; rigorously so, because the momentum-dependent terms in the 
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Figure 1: Notation for one- and two-gluon vertices. 



fermion propagato! are bounded. The behavior is the same for a closed fermion 
loop with any number of gluons attached. Hence, internai heavy-fermion loops 
decouple precisely as expected. 

The self-energy and vertex corrections are less trivial, because the external 
momenta are set on shell. Fig. Q shows the one- and two-gluon vertices. For the 
action S = So + Sb + Se (with the four-leaf dover for Se) 

Aq^ = 70 cos(p + ^/c)oa — i sm{p + ^k)oa 

(8.3) 

+\cEC,(yoj cos iA;oasin/!;ja, 
= Cli cos(p + \k)ia - ir sCsìi\{p + \k)ia 

(8.4) 

— ^CEQcTQj COS \kja sin k^a + ^CBC^^ij cos ^/cja sin kja, 

where p and k are the incoming fermion and gluon momenta, respectively. The 
expression for X^^ is not needed, except to note that its mass dependence is 

qualitatively the same as aI?^ . 

A formai way of going to the mass shell is to put po = ìE^q, with E from 
eq. ( [4. 41) . The 70 in the analytic continuation is not rigorous when applied under 
an integrai, but the mga dependence comes out right. The temporal vertex aH^' 

is proportional to e^' = 1 + moa. The spatial vertex Af^^ is proportional to Cj 
which, when it is tuned so that Mi = M2, satisfies r^C ~ (1 + moa)/Mia for 
large mass. Similar behavior holds for quark-multi-gluon vertices. Finally, the 
inverse propagator is also proportional to 1 + moa, for p dose to the mass shell. 

Consider any process with an external fermion line. Loop diagrams can be 
built up from the tree diagram by adding more gluons. Each additional vertex 
on the external line requires an additional fermion propagator. The dominant 
mass dependence of the propagator-and- vertex combination is (l + mo)/(l-|-mo) 
or C/(l + "T-o)) ^-nd thus cancels always. 
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For example, ali diagrams in the self-energy are proportional to 1 + itiqU. 
After summing the geometrie series and integrating over po one finds 



l + 5gM™(moa) + ••• 



[moa) 



Mìa 



1 



(8.6) 



where M\^\mQa) and ^^^^(moa) depend mildly on the mass, varying smoothly 
from the vaine obtained for massless fermions to the vaine in the static formn- 
lation. 

The same happens to the fermion-glnon vertex. The gauge-conpling renor- 
malization factor is defined throngh the fermion-glnon vertex via 



ZgJ\f{p)ux^t{p)K^{p,p)ui^t{p)J^{p) = -^U^eì{phf,Ureì{p) 



E 



•7) 



where A^{p',p) is the full vertex function, including leg contribntionsj^ In 
pertnrbation theory one usnally organizes the calcnlation by treating the legs 
and the proper vertex separately. By gange invariance 



Z, 



ZlF 



where Z3 (Z2) and Zi (Zip) are the gluon (fermion) wavefunction and proper 
vertex renormalization factors. The strong mass dependence of Z2 must, there- 
fore, cancel against Zip. (The residuai mass dependence of Zif/^2 should be 
the same as Z1/Z3 to satisfy the expectations of decoupling.) Indeed, at tree 
level the temperai vertex provides the asymptotic factor 1 + moa, and, by the 
general argument, the full proper vertex has the same (dominant) mass behavior 
to ali orders. Hence, 



Z^p{moa) 



1 - 9oZ^^j{moa) + ... 



(8.9) 



where z|p(moa) again depends only mildly on the mass. With the spatial ver- 
tex aI"' the factor compensates for the missing factor of 1 -|- nioa to ensure 
that I/M2 appears, so the 1 -|- niQa counting is the same. 

For electroweak currents and four-quark operators, the analysis of the mass 
dependence is similar. Again, loop diagrams have the same leading mass de- 
pendence as tree diagrams for the same process. For example, the bilinear Jp^, 
defined in eq. has renormalization Constant 



Zr{mfa, mgo) = e 



(Mija+A/iga)/2 



1 + glZ^^^ {mofa, mogo) + 



(8.10) 



^^At tree level one verifies = 1 from AP' and also from A''^' if ma — M2 
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The mass dependence of the loop corrections (mo/a, niQga) smoothly connects 
massless and static results. Such behavior is borne out in sect. |l^'s nonpertur- 
bative check of the locai vector current, for which F = 79. 

The considerations of this subsection argue that the large-mass limit of ac- 
tions described by eq. is well-behaved in perturbation theory. More gener- 
ally, the physical masses and, hence, the couplings could depend on the gauge 
coupling in a nonperturbative way. But because the origin of the gauge-coupling 
dependence is the region of momentum space near the cutoff, it seems unlikely 
that nonperturbative contributions would overwhelm the perturbative contri- 
bution, at least once the cutoff is large enough. Should perturbation theory 
prove inadequate, however, a nonperturbative renormalization group could, in 
principle, substitute for perturbative calculationsj^ Nevertheless, it seems im- 
plausible that nonperturbative effects are more worrisome at large mass than at 
small. Thus, the main conclusion, that the large-mass behavior of interacting 
fermions is benign, is probably valid nonperturbatively. 

8.2 Mean field theory 

To estimate the one-loop corrections, recali that the dominant contributions 
come (in Feynman gauge) from tadpole diagrams, which originate from higher- 
order terms in the expansion of the link matrix C/^ = 1 + go^/x + ^9o^fi + ' ' It 
is possible to make this observation more systematic [ ^] . Wherever the gauge 
field appears, substitute 

U^ix)^uo[U4x)/uo], (8.11) 

where uq is a gauge-invariant average of the link matrices. The substitution 
should be understood in the following sense: The second factor [Ufj,/uo] produces 
perturbative series with small coefRcients. The first factor uq, which has a nasty 
tadpole-dominated perturbative series, should be absorbed info the couplings c„ 
and into renormalization factors Zo. When a numerical value for uq is needed, 
for example in a Monte Carlo calculation, it should be taken from the Monte 
Carlo itself. 

With this prescription the hopping-parameter form of Sq remains as in 
eq. (|1.6|), but with Uf^ U^/uq and 

K-s^t l^s,t = UQKs^f (8.12) 

The mass form of is given by eqs. (|2.3| ) but with difference operators defined 
with 

f±^ = no ^T±^, (8.13) 

instead of T±^, and mass 
^^For example, to tune nonperturbatively, compute the energy of a meson and imposing 

Mq = M2. 
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moa = _ [1 + r,C{d - 1)] = !^ + [1 + r,C(d - 1)] U,' - l) (8.14) 

instead of mo- Finally, an overall factor of uq multiplies each term in the action. 

The clover-leaf construction used to define the chromomagnetic and chromo- 
electric fields contains products of four U matrices. If one replaces the gauge 
fields B and E with tadpole-improved clovers, the interactions Sb and Se are 
given by eqs. (|1.7| ) and respectively, but with 

CB = U^CB, CE = U^CE, (8.15) 

instead of and ce, and kg instead of Kg. The fourth factor of uq corresponds 
to the overah factor mentioned above. 

After these rearrangements one can immediately generahze the expressions 
in sects. I and I to the mean-field level. They remain the same as before, 
but with niQ — > mo, cb — > cb, and ce — > ce- Consequently, the couphngs 
cb, and ce{0) should be adjusted to the right-hand sides of eqs. ( |5.32| )-( 5.34| ), 



but with mo ffiQ. The resulting conditions represent a set of mean-field- 
theory predictions at Qq ^ 0, given a nonperturbative input for uq. One-loop 
calculations with mo 7^ 0, 7^ 0, and c^; 7^ wih test and correct mean-field 
theory estimates. 

At currently accessible lattice spacings ref. [ has shown that, with this 
mean-field reorganization and a sensible choice of expansion parameter, the 
bare perturbative series converges quickly in many cases. Calculations [ 18| 



in one-loop perturbation theory of Feynman diagrams needed to determine the 



[1] 



eh 



(moa) show a smooth transition from the massless to the static limits.PI One 



therefore expects the essential concepts of ref. [ to apply to the Cn and to the 
coefficients in eq. ( [7.1| ) too. Indeed, in the one case for which a nonperturbative 
check is unambiguous, the normalization of the vector current, there is excellent 
agreement with mean-field theory, cf. fig. in sect. IC. 



9 The Nonrelativistic Limit 

It is illuminating to adapt the methods of sect. ^ to the nonrelativistic and static 
limits. Rather than adjusting the couplings to obtain the Dirac Hamiltonian, 
one could instead aim for the nonrelativistic Pauli Hamiltonian (and generaliza- 
tions thereof). An advantage of this avenue is that it provides a useful physical 
picture even when the couplings are maladjusted, in particular when ^ = 1. 
Many Monte Carlo studies have used actions with C = li and it would be help- 
ful to have a framework for interpreting their data in the heavy-quark regime. 
Indeed, the analysis of this section shows that C = 1 is acceptable for nonrel- 
ativistic fermions. Even the Wilson action (somewhat crudely) approximates 

^^For tadpole and scale-choice improvement [ ^ of the static limit and of nonrelativistic 
QCD, see refs. ||. 
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the properties of nonrelativistic or heavy-quark effective theory, provided niQa 
is adjusted correctly. Similarly, the Sheikholeslami-Wohlert action is a better 
approximation. 

The Hamiltonian of the action 5o + Sb + Se can be brought to the nonrel- 
ativistic Pauh form with the Foldy-Wouthuysen-Tani transformation. Imagine 
transforming the -f D and a- E terms away completely. Afterwards, the trans- 
formed Hamiltonian reads 

where Mi, M2, and Mb are as in sects. ^ and ^. The new mass Me reduces 
to M2 with suitable mass dependence of ce (cf. Appendix or as rn-oa 0. 
The specific expression is not needed here. The Pauli form of eq. ( |9.lD has 
no coupling between the upper (particle) and lower (anti-particle) components 
of as in the explicitly nonrelativistic formulations [ |l2|. Here, however, 
eq. ( |9.1| ) is derived within the lattice theory, rather than being an Ansatz for an 
effective lattice theory. 

Let US discuss the physics of each term in eq. (|9.lD . The first three are the 
rest mass, Coulomb potential, and kinetic energyfn of the fermion. The S • B 
term, as one recalls from atomic physics, produces the hyperfine splitting. The 
last term can be rewritten 

[■y-D,j-E] = iT,-{Dx E-Ex D) + {DE - ED). (9.2) 

The two parentheses give the (non-Abelian) spin-orbit and Darwin interactions, 
respectively. 

The Pauli Hamiltonian is quantitatively useful only if the fermion is nonrel- 
ativistic. Given nonrelativistic velocities, however, eq. ( |9.1| ) remains applicable 
even when the various masses are unequal. Fig. |2| is a sketch of the quarkonium 
spectrum, illustrating how the masses affect the spectrum. The interesting gross 
feature of the spectrum is not the overall mass gap — dose to 2Mi — but the pat- 
tern of radiai and orbitai excitations, e.g. m2s — "iis or mip — mis. These 
splittings are dictated by the kinetic mass M2. Following the analysis of ref. [ 
|l2|| they are of order M2f^, where v is the typical velocity of a heavy quark in 
quarkonium. {v 0.3 for charmonium, and v ~ 0.1 for bottomonium.) Further 
application of the velocity counting in ref. [ |l2| to eq. (|9.1| ) shows that the hy- 
perfine splittings are H-B/Mb ~ Mlf^/M^, and the spin-orbit splittings are 
[7.D,7.£;]/M|~M2V/M|. 

The preceding paragraph merely reviews the well-known argument that the 
rest mass of a nonrelativistic particle decouples from the interesting dynamics. 
In our formalism the reasoning suggests the following strategy: forget about Mi 
and adjust the bare mass so that the kinetic mass M2 takes the physical value. 

^"^Because of this physical interpretation the quantity M2, defined in eq. ( |4.7| ), is called the 
kinetic mass. 
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Figure 2: Quarkonium spectrum and the influence of the masses Mi, M2, M^, 
and M^;. (A similar picture apphes to the heavy-hght spectrum, except the 
overall gap is Mi instead of 2Mi, and orbitai and radiai excitations are set by 
Aqcd instead of M2V^ .) 



Meanwhile, choose the coupling C, by convenience. The obvious example is to 
take C = 1, as in the Wilson and Sheikholeslami-Wohlert actions. 

Since the Wilson and Sheikholeslami-Wohlert actions represent viable non- 
relativistic field theories, it makes sense to compare them to the explicitly non- 
relativistic theories. The (tree-level) masses for the Wilson action are plotted as 
a function of rriQa in fig. ^. Assuming rriQa is chosen so that M2 = f^q, the other 
masses satisfy Mi < rug, M^^ < , and M^^ < m~^. The simplest form of 

nonrelativistic QCD [ 11| has Hamiltonian Hj^n = ^D'^'^ /{2mq). Thus, in our 



notation, -ffNR has Mi = and M^ = M^ = 0. Thus, the Hamiltonians of 
the Wilson and simplest nonrelativistic theories make the same errors qualita- 
tively. For example, in both one expects the fine and hyperfine splittings to be 
too small. Similar ly, for the Sheikholeslami-Wohlert action one finds Mb = M2, 
and thus good hyperfine splittings, but M^^ > m~^, so the splittings between 
Xj states ought to be too large. To obtain the correct spin-orbit splittings, one 
needs the mass dependence of c^;, cf. Appendix 

A parallel set of remarks applies to heavy-light systems. The Hamiltonian 
of the lattice theory satisfies the usuai heavy-quark symmetries as niq — > cxd, 
no matter what M2, M^, and Me are. On the other hand, the lattice theory 
possesses the rightP| l/mq corrections only if M2 = Mb = mq. A computation 
with the Wilson action and mq = M2 obtains spin-averaged features correctly, 
but underestimates the chromomagnetic l/w-g corrections. Compared to cor- 

^^Many phenomenological applications require matrix elements of operators of the elec- 
troweak Hamiltonian. These operators must also be constructed to the appropriate order 
in 1/mQ, cf. sect . | ^ In particular, to first order in I/mg the coefRcient di must be chosen 
according to eq. (7.1C). 
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Figure 3: The (tree-level) masses for Wilson fermions (C = 1, Ts = 1). By 
happenstance M2 is always within 15% of niQ, which is a result of a conspiracy 
between the 'y-D and Wilson terms. 



rections from the kinetic energy, the spin dependent effects are thought to be 



small [ 17], so again the most essential adjustment is M2 = mQ. A better com- 



putation with the Sheikholeslami-Wohlert action and mg = M2 = Mb obtains 
all^^ l/mq features correctly. 

Despite the similarity between previous nonrelativistic field theories [ ^, ^, |l^, 



11, 12| and the view adopted in this section, there are significant technical differ- 



ences. The four-component approach explicitly includes the terms mQ^lJ(x)^lJ{x) 
and terms that couple upper and lower components, such as iJj{x)-f-Dip{x) 



and 'il){x)ciL-Eij){x). The program of Lepage, et al, [ 12| omits these interac- 
tions in practice, though perhaps not in principle. There are advantages to 
leaving out the Dirac block-off-diagonal interactions. Fermion propagators are 
the solution of a (one-sweep) initial-value problem, whereas they are otherwise 
the solution of a boundary-value problem, solvable only by iteration; with fewer 
interactions, perturbation theory is easier [ |2^. On the other hand, these in- 
teractions are necessary to take a ^ (by brute force) without the scourge of 
power-law divergences, or to reach info the semirelativistic regime. 
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10 Numerical Tests and Applications 



With a few examples this section tests the results of the previous sections with 
Monte Carlo data. Ali data were generated with the axis-interchange symmetric 
Wilson or Sheikholeslami-Wohlert actions, so as m^a increases, we rely on the 
nonrelativistic interpretation of sect. ^. The tests verify the most important 
lessons. The bare mass should be adjusted until the kinetic mass M2, defined in 
eq. ( [4. 71) , takes the desired vaine. In particular, in an extrapolation to vanishing 
lattice spacing, one ought to hold the kinetic massf^ fixed. On the other hand, 
the dynamically irrelevant rest mass may deviate from mg. For matrix elements, 
it is also important to use the improved field of eq. (7.8) or ( [A. 17 ). The factor 



gMia/2 ]^Qj,g important than the bracket in eq. ( [7.8D , because it guarantees a 
smooth approach to the static limit.|^ 

First, consider the mass spectrum, in particular the hyperfine splitting in 
heavy-light systems. By heavy-quark spin symmetry the vector-pseudoscalar 
mass difference my — mp is expected to be proportional to I/mg. Obviously 
the leading term in the sum my + mp is proportional to mg, so the combination 



ruy — rup should be nearly independent of mg. Numerical work [ 29 1 found, 
however, that rriy — uip decreases for increasing mg, with lattice spacing a 
held fixed. These analyses take my and mp from the rest mass. From sect. ^ 
the rest mass Mi of the quark governs the rest mass of the mesons, while the 
chromomagnetic mass Mp governs the hyperfine splitting. The computed lattice 
quantity, therefore, is proportional to Mi/Mp, which decreases for increasing 
quark mass, cf. fig. ^. Given Mi, Mp is not as large with the Sheikholeslami- 
Wohlert action as with the Wilson action. Numerical data with cp = 1 show 
behavior [ qualitatively similar to cp = 0. 

To improve the determination of rriv — mp one should tune to the kinetic 
mass instead of the rest mass and use the mean- field or one-loop estimate of cp. 
The chosen vaine of cp could be tested in quarkonia. Then in heavy-light systems 
one could verify two predictions of heavy-quark symmetry, as applied to the 
lattice theory: a falling Mi/Mp behavior when using the mesons' rest masses, 
and a fiat M2/MP behavior when using the mesons' kinetic masses. 

Next, consider the improvement and normalization of multi-quark operators 
from sect. 0. The normalization of the vector current can be checked nonper- 
turbatively. The fermion number 



($P$T^) 



counts the number of /i-flavored fermions in \H). If \H) has one and only one h 
in it, the condition Nh = 1 defines the factor 2KtZv^ Fig. ^ compares this non- 

^^This can be done nonperturbatively with a meson instead of a quark state. 
^'^Neglecting the bracket introduces only 0{mopa^) lattice artifacts at moa <C 1, but 
0{p/mq) at moa ^ 1. 

^®In the notation of sect. tì, Zy = Z^^ and Za ~ Zy^^^. 



37 



0.8 - 



0.6 - 



0.4 



0.2 



P = 5.7, c = 0.0, Kc = 0.16913 
|3 = 5.9,c = 0.0,Kc = 0.15974 
|3 = 5.9,c= 1.4, Ke = 0.14159 
P = 6.1, c = 0.0, Kc = 0.15495 



T □ 



▼ □ V 



0.2 



0.4 



0.6 



0.8 



mo 



1 + mo 



Figure 4: The charge normalization factor 2kZv vs. (kc — k)/{hìc — 3k/4) = 
(1 — 8k)/(1 — 6k) = m-o/(l + ?7T,o), with uq = 1/8kc. The symbols are Monte 
Carlo determinations with = C = 1 ^-iid = c^; = c. The sohd square is 
the exact result for k = 0. The sohd curve is the mean-field approximation to 
eq. ( |7.7D , 2K,Zy = 1 — 6k. The dashed curve is a mean-field Ansatz 2R, which 
(foolishly) neglects the mass dependence. 



perturbative definition of 2kZv with the mean-field-improved, tree-levei pertur- 



bative approximation. The symbols are from Monte Carlo calculations [ 3C] of 
eq. ( lO.ip , with \H) a meson with a spectator anti-quark of different flavor, and 
the solid curve is the mean-field improved, tree-level approximation 2K(l + moft,)- 
Fig. ^ exhibits several interesting features. The solid curve accurately tracks the 
dominant mass dependence from friQ = to mo = co. From eq. ( ^.10| ) one ex- 
pects a subdominant mass dependence from loop corrections Zy {rrioha). Indeed, 
near mg = the massless one- loop correction [ ^] accounts quantitatively for 
the discrepancy, and near mo = oo the discrepancy becomes smaller, in accord 
with a Ward identity, which requires 2kZv = 1 at infinite mass [ |3^ . Neglecting 
the dominant mass dependence, as in the dashed curve, is obviously completely 
wrong for mg ^ 1 . 

Finally, consider the decay Constant of a heavy-light pseudoscalar meson, 
computed with the locai axial current J^^{x) = '(/'"(x)7^75V'''(x). Fig. ^ shows 
Monte Carlo data [ |3|, || at /? = 5.7 (for which a'^ ^ 1 GeV) for 
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Figure 5: (j)p = ^mpfp vs. 1/mq. Squares represent resulta from the con- 
ventional zero-mass normalization, plotted versus 1/Mi. Circles represent the 
same results with the correct normahzation, plotted versus I/M2. The points 
for > are obtained with the Wilson action [34]. The static {m~^ = 0) 
point is from Ref. [ E5|. The curves guide the eye, and the approximate location 



of the physical B and D mesons is shown. 



cl>p = V2 {0\Zfj;il\P, 0) = ^fpVf,, (10.2) 

where is the meson's four-velocity. (The vacuum and one-meson states are 
normalized to unity.) We have deliberately chosen a largish lattice spacing to 
enhance lattice artifacts, and thus test our control over them. We analyze the 
data two different ways. The lower set of points takes the meson mass from the 
rest mass and neglects the factor-^^ Z"^ = (|7.8D . As suggested 

by the curve, the neglectful analysis would produce a locus of points that ap- 
proaches zero in the static limit. The upper set of points uses the normalization 
factor, and — ^just as important — it defines the meson mass through a mean-field 
approximation to the kinetic massj^ Fig. ^ shows how cruciai both refinements 
are, if the Wilson-action data are to approach the static limit smoothly. 

An important application of a plot like fig. |5| is to compute the slope in 
1/mq. In the heavy-quark effective theory the slope is seen to arise from three 
sources: the kinetic energy, the chromomagnetic interaction, and a correction to 

^^Ref. [ Q provides hopping parameters and rest masses only, but the mean-field approxi- 
mation is adequate for illustrative purposes. 
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the infinite-mass current 
energy requires tuning M2 
tuning Mb = ma, 
compactly by eqs. 
the correct slope [ 



17]. The lattice theory has direct analoga: the kinetic 
the chromomagnetic contribution requires 



11 



and the locai lattice current requires a correction, given most 
I) and ( 7.10|) . Ali three ingredients are needed to obtain 



11 Conclusions 



This paper (and conference reports [ 3C] anticipating it) provides the 

foundations of a theory of lattice fermions, valid at any mass — large or small. 
The action starts with a set of interactions that encompasses those of both light- 
fermion actions [ 13, ^ and heavy-fermion actions [ 10, 11, 12]. The couplings of 
such a general action are then tuned in successive approximations to the renor- 
malized trajectory. In applying renormalization-group techniques to analyze and 
reduce cutoff effects we do not, however, expand in either m^a or Agco/^ig. 

Although there are several methods for tuning the action (for example, refs. [ 

|39|]), our analysis is based on Symanzik-like, on-shell improvement criteria. 
This entails the computation of on-shell correlation functions, or, equivalently, 
of the Hamiltonian. Enforcing continuum-limit behavior — for example the rela- 
tivistic mass shell — yields conditions on couplings of higher-dimension interac- 
tions. In practice, we here compute on-shell quantities in (tadpole-improved) 
perturbation theory. 

An examination of the lattice theory's Hamiltonian, derived from the trans- 
fer matrix, is especially illuminating. It shows that it is unnecessary to im- 
prove Wilson's discretization of the time derivative. Instead higher-dimension 
interactions can be built from fields on one (or, in some cases, two adjacent) 
timeslices. Thus, our class of improved lattice actions automatically has an 
easy-to-construct transfer matrix. The actions, consequently, ali automatically 
satisfy heavy-quark symmetries in the limit of large mass. 

The Hamiltonian is a useful tool for examining lattice artifacts. A term of 
dimension s + 4 factorizes as 

5Hut = a'b{mqa,g^)H, (11.1) 

with dependence on the theory's relevant couplings, m^a and g^, in the coeffi- 
cient b. Prom this expression it is plain how the artifacts behave as the mass 
increases: the absolute error induced by eq. ( |11.1 ) is {5H\^i) ~ {paYp, where p 



is the typical three-momentum of the process. At small mass it is standard that 
the associated coefìicient is a benign numerical factor 6(0). At large mass the 
lattice action's heavy-quark symmetry implies that it is a (generally different) 
numerical factor b{oo). Our explicit results at tree level and our analysis of 
higher orders show that b{mqa) is a smooth, gentle function connecting the two 
extremes. In a nutshell, therefore, the characteristic measures of cutoff artifacts 
are Aqcdo and pa, but never m^a. 
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In general our actions have two hopping parameters. Then it is possible to 
maintain equality between the rest mass (energy at zero momentum, eq. ([4.6|)) 
and the kinetic mass (inertial response, eq. ([4.7|)). In nonrelativistic systems, 
however, there is a noteworthy simphfication. Embracing the philosophy of the 
static [ ^ or nonrelativistic [ |ll|] theories, one can ignore the rest mass 
and, hence, forgo one of the hopping parameters. The obvious application is to 
set them equal, as in the Wilson [ ^ and Sheikholeslami-Wohlert [ ^ actions. 
Therefore, the correct interpretation of numerical data generated with these 
actions at m^a ~ 1 is a nonrelativistic one. In particular, the hopping parameter 
must be adjusted so that the kinetic mass agrees with the physical mass. 

If one knows a priori that a quark is nonrelativistic, it is computationally 



cheaper to use a two-component formalism [ 12|. But there are many instances in 
which one would like to trace the mass dependence from the static limit down to, 
say, the strange quark. One example is fig. |5|, which, with reliable calculations, 
should indicate how and where the heavy-quark expansion deteriorates. 

The results of this paper can be extended in several ways. The couplings have 
been computed at tree level with mean-field improvement. One-loop calculations 
are desirable, and better stili would be a nonperturbative determination, perhaps 
in a mass-dependent generalization of ref. [ 40 1. Once one is confident that 
0(a) artifacts are under control, one can extend the analysis to dimension-six 
interactions. Tree-level, O(a^) improvement should be manageable; beyond tree 
level the bothersome four-fermion interactions enter the fray. 

One would like to use the actions presented here in Monte Carlo calcula- 
tions of QCD. If one uses the 0(t;^)-improved action to compute the spectrum 
of charmonium and bottomonium, then, without more tuning, one could calcu- 
late properties of D and B mesons, including the electroweak matrix elements 
needed to determine the unknown elements of the Cabibbo-Kobayashi-Maskawa 
matrix.0 Note that for these matrix elements, as well as for the quark mass 
in the MS scheme, a small lattice spacing is helpful to reduce perturbative cor- 
rections. Once the lattice spacing is small enough so that mga^l, our for- 
mulation is especially advantageous. The two-component nonrelativistic theory 
breaks down as mqa gets smaller [11], yet the old small-mass theory would have 
leading lattice artifacts of order asmqa and (mgo)^. Our improved action, on 
the other hand, remains viable for any mass, and its cutoff effects are small, of 
order a^AgcDO and (Aqcdo)^- To obtain comparable accuracy through brute 
force in the old theory, one would have to reduce the lattice spacing by a factor 
of mg /Aqcd — about live for the charm quark. Even for a perfect algorithm the 
savings in computer time is, therefore, a factor of 5*^. 

20pQj. j-elevant reviews, see ref. [ ^ . 
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Table 5: Estimates of the size of the dimension-seven mteractions that arise 
in designing an action for quarkonia with an accuracy of mqv^. (There are 
many other interactions needed to ensure ©((Agcoa)^) accuracy in all-hght and 
heavy-hght systems.) 



Hn 


only hght 


heavy-hght 


quarkonia 


^(1)2)2^ 


^^QCD 


A4 

^^QCD 


4 4 

niQV 




A4 

^^QCD 


A4 

^^QCD 


4 4 
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A Quarkonium to 0{v^) 

This appendix extends the analysis of the main text to incorporate interactions 
that contribute in quarkonium through O(f^). Naively, this would entail dose 
scrutiny of ali interactions through O(p^), i.e. up to dimension seven. Some 
dimension-seven bilinear interactions are listed in Table ^, with their magnitude 
in quarkonium estimated by the velocity-counting rules of ref. [12]. Together 
with Table g one sees that not ali dimension-six and -seven interactions are 
necessary to Oiv'^); this Appendix considers only the entries that are. 
One must consider the dimension-six interaction 

Sso = CsoHs^Ì)nlo[l-D,l-E]ll)n, (A.l) 

n 

with coupling Cso, and the dimension-six and -seven interactions 

^4 = -2ks [c4Dl-D - ir,C4LA(3)] A^^)^.^, (A.2) 

n 

with couplings C4D and 04/,. Further dimension-six and -seven interactions con- 
tribute in 0(^;^) or higher [12]. We discuss the adjustment of Cso in eq. ( [A.l ) 



in sect. |A.1| , and the adjustment of C4D and C4L in eq. ( [A. 21) in sect. A.2| . 



The discretization of the covariant difference D and Laplacian A^^^ also must 
be improved, to remove ^^iDf"^ and ^Df"^, respectively. These interactions 
break rotational invariance, and we treat them in sect. |A.2| . 
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A.l Chromoelectric interactions to O(f^) 

It is easiest to treat the "spin-orbit" interaction 5*80 in the Hamiltonian formahsm 
of sect. ^. The coupUng Cso can only appear in a dimension-six terni imphed by 
the elhpsis in eq. ( 5.16D . The spin-orbit Hamiltonian is 

Hso = a\o{moa)4>jo[j-D,'yE]^. (A.3) 
Under the two-parameter Foldy-Wouthuysen-Tani transformation of eq. ( p.23| ) 
Ko = bso - èCi + bE^i + biÌE - 2mgabo^iCE- (A.4) 
With òso it is possible to give an invariant involving Òe- 
Be = bi—4mqabQbibE—8{mqabQ)'^bso = b'^ —Amqab'Qb'ib'^ — S{rnqab'Qf'b'g„, (A. 5) 

and one wants Be = 1- Just as the redundancy associated with the Foldy- 
Wouthuysen-Tani parameter intertwines the mass dependence of C and r^, 
the redundancy associated with the other Foldy-Wouthuysen-Tani parameter S^e 
intertwines the mass dependence of ce and Cgo too. 

Assuming a discretization of 5so that resides on two timeslices only, it is 
straightforward to generalize the transfer-matrix construction to the action So + 
Sb + Se + Sso- After expanding the transfer matrix in powers of a one finds 

òso = - CE)ef2{mo) + ^TT^T^- 

Z[i + mo) 

Combining eqs. ([A^), (|]2|), and (|AJ) and setting Be = l yields the mass- 
dependent condition 

^ - 1 rsC r^mo(2 + mo) Csomo(2 + mg) 

'^^ ~ mo{2 + mo) 1 + mo 4(l+mo)2 C(l + "^o) ' 

Here the rest mass Mi has been eliminated in favor of the kinetic mass M2 as 
appropriate to the nonrelativistic interpretation of sect. ^. 

The redundant direction associated to S,e permits a free choice of Csq. In our 
framework, which stresses a smooth matching to the massless limit, the most 
convenient choice is probably Cgo = 0. But for purely nonrelativistic applications 
ref. [ ^] would choose ce = and Ccso oc m~^. Other possibilities correspond 
to the special choice for in eq. (|5.35|) . Then eq. ( A.7 ) reduces to 

CE = l + 21og(l + (A.8) 

The further special case corresponding to = is c^; = 1 (independent of mga) 
and Cso = 0. 



43 



A. 2 Kinetic energy to 0{v^) 

The interaction S/j^ produces corrections to the kinetic energy. It is easiest to 
analyze from the energy-momentum relation, as in sect. Expanding eq. ( |4.4[ ) 
to O(p^) yields 



where 



and 



/ a4p \-V3 



1 d^E 



4 dp: 



n 4M| ■ 

p=0 * 



(A.ll) 



The relativistic mass sheh satisfies M4 = M2 = Mi, and a nonrelativistic mass 
shell with leading relativistic correction satisfies M4 = M2. In both cases rota- 
tional invar lance requires 11)4 = 0. 

A straightforward way is enforce t/;4 = is to take an improved covariant 
difference 

aA = im - T-^) -UTf- Tìù (A.12) 
and an improved covariant Laplacian 

^2^(3) ^ ^ U^rj..^ ^ ^2)~^{T^ + . - 2)) . (A.13) 



The coefficients are chosen so that the Fourier transforms have no or ^ipf 
terms, respectively. Then (at tree-level) u;4 = automatically. 

Since S4 contains no higher time derivatives, the transfer-matrix construction 
proceeds as usuai. After deriving M4 for Sq + 5*4 (Mi and M2 are unchanged), 
one finds M4 = M2 at tree level if the couplings c^d and c^l obey 

4C2c4£,(l + mo) + rsCc4Lmo{2 + mo) = 

C^(l + mo)[2(l-C^) + mo(2 + mo)] r,C^[2(l + mo)^ - 3C^] 

m2(2 + mo)2 + mo(2 + mo) ^^'^^^ 

rgC^[mo(2 + mo) - 6C^] _ r^C^mo(2 + mo) 
4(1 + mo) 4(l + mo)2 

This result holds whether is tuned so that Mi = M2 or not. As with C and 
Ce, only the massless limit of 04/5 is unambiguous. For the full mass dependence 
the redundant 04^ must be specified, for example c^l = 0. 
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It is instructive to look explicitly at the consequences of omitting 5*4 from 
Monte Carlo calculations. With wnimproved D and A^^^ and no 5*4 



and 



- mo(2 + mo) + 4(1 + mo) ' ^^'^^^ 



8C^ ^4C3[C + 2..(l + ^o)]^ ..V ^^^^^^ 



M| m|](2 + mo)3 mg(2 + mo)2 (1 + mo)^' 

The rotational-invariance breaking artifact is, thus, Oipfo? /mq) for m^a large 
and small. The rotationally invariant, relativity-breaking artifact is 0(j)'^a? /niq) 
at small m^a, and ©(p^a^/m^) at large m^a. 

A. 3 Electroweak operators 

For electroweak decays of quarkonia to O(v^), one needs a higher-dimensional 
generalization of eq. (^), 

^^[x) = e^i"/2 h _^ a^^^^.i:) + ia2d2A(3) 

(A.17) 

+ ^a^dBT.-B + \a'^dEOL-E 

The d's are easiest to derive from the Foldy-Wouthuysen-Tani transformed field. 
Combining eqs. dJ) and (|]23|) 

■^{x) = exp(aa7-D + a'^ÌEOi-E)e^/'^i;{x), (A.18) 

where and parameterize the solution of the tuning conditions. Expanding 
the cumbersome exponentials to 0{a?'\p'^ , B , E]) and eliminating and ^e iii 
favor of the couplings C-, i~s, cb, and ce, one finds 



d2 = di 



2 ^sC 



2(1 + mo) 

cbC 
2(1 + mo) 
C(l - ce){1 + mp) _ 
mo(2 + mo) M2 ' 



2(1 + mo) 

r^^ - 

dE ■■ 



and di as in eq. (ItIcI ). In eqs. ( |A.19D the kinetic mass M2 has been substituted 
for the rest mass Mi. Thus, these formulae remain valid under the nonrelativistic 
interpretation explained in sect. 0. 
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B Combining Exponents 



This appendix presenta a proof of eq. ( [5. 261) , i.e. that the functions fi and /2 
are given by the expressions in eq. ( 5.28| ). If Hq, Hj, and Hj were to commute, 



it would be trivial to combine the exponents. They do not, so the combined 
exponent depends on their commutators as weh. The commutators are 

[Ho,Hi] = -2MiHj + 0{p^a^), 

[Ho, È}] = 2MiH\ + 0(p3a3), (B.l) 

[Hi, h\] = -C^l-G^^r + 0(p4a4). 

To obtain the last commutator we have written Djj = [D — ^E)cont and 
Dy = {D + i£!)cont5 and we have neglected higher powers of -Econt- The op- 
erators Q and A^^^ carry one and two powers of pa, respectively; thus M. = 
Mi - ir,Ce-^iA(3) +0(p4o4). Note that although [Hi,h\] is Oip^a^), further 
commutators such as [Hq, [Hi,h\]] are at least O(p^a^). 

Instead of solving the field theory, it is enough to consider a toy model 
with two degrees of freedom, a fermion (annihilated by a) and an anti-fermion 
(annihilated by h). With discrete time the action is 

S = Y,a\ (do + m)at + b\ {8^ + m)bt - ii9 (a\b\ - hat) . (B.2) 



The transfer matrix has the same form as eq. ( p.25 ) with 



Ho = log(l + m) (iU + È^È^ , Hi = i^ÈÀ, (B.3) 

where A = {1 + m^/'^a, and B = {1+ rriY^'^b. With the identification of -d with 
C,Q and m with e-^ — 1, the operators and È'f^ of the toy model and the 
field theory have the same algebraic structure. 

This model has only four states, the vacuum, fermion, anti-fermion, and a 
fermion-anti-fermion state. The strategy is to work out the transfer matrix 
elements explicitly, and then take the logarithm. These steps are easier in the 
Grassman-number approach, where the matrix elements of the transfer matrix 
are the coefficients of monomials in (Grassman numbers) A, A^ , B, and B\ 
when 

T{AIB^;A, B) = ei'5«tbt^ata+6t6g-itìba 4) 

is expressed as a polynomial. Up to factors analogous to det(2Kji?), which we 
can drop without loss, the transfer matrix in the neutral sector is 
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Writing T = VDV\ where D is diagonal, the Hamiltonian is H = —V ìog{D)V^ . 
Expanding the result to ©("ì?^) one finds 

where e™^ = 1 + m. Expressed in terms of Fock-space operators 

H =[mi- /2(m)i?2] i^^ti _ ^^t^ _ if^irn)^ (À^ - Si) . (B.7) 



Substituting mi and d for M. and completes the derivation of eq. (5.26). 



C Spinors, Creation and Annihilation Operators 

This Appendix gives the construction of spinors and of creation and annihila- 
tion operators in d = 4 space-time dimensions. These are needed to calculate 
amphtudes of on-shell fermions via Feynman diagrams. 
Consider an arbitrary bihnear fermion action 

5 = E V^(^) (l^^K^ix, y) + Z(x, y)) V'(?/) (CI) 

with an imphed sum over /i. We assume that and L are translation in- 
variant. With parity L{x,y) is symmetric, and Kf^{x,y) anti-symmetric, under 
interchange of x and y. The field ^p{x) has the following equation of motion 

E [j^K^ix, y) + L{x, y)) ^{y) = 0. (C2) 
y 

In such a free theory, one searches for solutions of the form ip{x) = e*^°*^*'''*n(p). 
The four-component spinor u{p) must satisfy 

{i^^K^{p) + L{p))u{p) = ^. (C3) 

The Fourier transforms K^[p) and L{p) are real functions of p by parity and 
translation invariance. Multiplying by —ij^K^{p) + L{p) one sees that solutions 
exist only if 

K\p) + L\p) = 0. (C4) 

This mass shell coincides with the one derived from the propagator, as in sect. |^. 

The actions that we consider ali have the Wilson time derivative. Writing 
p = {po,p), one then has 

Ko{p) = sin Po, 

K,{p) = Ki{p) (C.5) 
L{p) = fj,{p) - cospo; 
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Ki{p) is an odd, and an even, function of p. Thus, solutions exist only 

when 

PO = ±iE{p), (C.6) 

where 



It is convenient to label the solutions of eq. ( p.3| ) by the sign in eq. ( p.6| ). For 
each sign there are two solutions, u{±l,p) and u{±2,p). Setting p = the 
equation of motion simplifies to 

sinh Mi (-70 sign ^ + 0) = 0, (C.8) 



where Mi = £'(0) = log[/i(0)]. Choosing 70 as in eq. ( |2.6D the four solutions at 
p = are 

ni(l, 0) = U2{2, 0) = mi-l, 0) = U4Ì-2, 0) = 1, (C.9) 

where the subscript is the Dirac index, and ali other components are zero. Direct 
substitution verifies that 

solves eq. (|C.3| ) for p / 0, if sinpo = ^sign^ sinhE' and L = fi{p) — cosh E{p). 
The denominator yields the normalization convention in eq.( p.l4| ). 

The two solutions with "negative energy" < 0) correspond to anti-particle 
states. As usuai we introduce 

v{^,p)=u{-C,-p), e =1,2. (C.ll) 

The spinors v obey the equation of motion 

i-ij^K^ip) + L{p))v{p) = 0, (C.12) 

which is solved by 

^ _^^^EER±±^y^^^O)^ (C.13) 

^ ^ v/2L(L + sinhE) ^ ^ ^ ' 

now with sin Po = +isinh£'. From now on we shall use u and v with ^ G {1, 2} 
and sin Po = +isinììE only. 

The spinors obey the conventional orthonormality properties 

uiC',p)u{C,p) = ,p)v{Ì,p) = 5^'^, 

(C.14) 

,p)v{Ì,p) = v{Ì' ,p)u{Ì,p) = 0, 
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where u = u'^jo ^^à. v = v^jq. Moreover, 

u{^',Phov{C, -P) = v{Ì', -p)lQu{Ì,P) = 0. 

In a relativistic theory E/m would appear here. 

The general solution to eq. (|C.2| ) is a linear superposition 



-i K'^'^) (C.16) 



d){Ì,p)v{Ì,p)e- 



■ipot—ip-x 



with sin Po = isinhi5^. The normalization factor N{p) is fixed below, after 
invoking this expansion for Hilbert-space operators. The operator-valued ex- 
pansion coefìicients and (i^(C)P) create particle and anti-particle states 

respectively: 

k(e,p)> = h\i,pm, m,p)) = d\i,pm, (c.i?) 

where |0) is the Fock state annihilated by ali b{^,p) and d{^,p). Assuming the 
vacuum is normalized to (0|0) = 1, the fermion states are normalized to 

{q{i\p')m,p)) = i27Tf6Ìp'-p)6^'^<PÌp), (C.18) 

and similarly for the anti-fermion state \qiS,,p)), if and only if the anti-commu- 
tator 

m',p'),bHtp)} = {27rf6{p'-p)6^'^<P{p), (C.19) 

and similarly for {d{£^' ,p'),d^^,p)}. 

The transfer-matrix construction provides the anti-commutation relation for 
ip{t,x) and ìp{t,x). Eq. becomes 



^{t,x) = J2j (03^^''"^""^V(P)'/V(t,y). (C.20) 

After inverting the Fourier series and evaluating the anti-commutators one finds 

^(P) -cosh.E 
•^(^^ = Kp)MHp)smhE- (^-''^ 

The convention (j) = 1 is the most convenient.p^ Then 

^^Another candidate is the relativistic convention, which is not at ali naturai in nonpertur- 
bative calculations. How can one normalize to = m? + /m without solving the theory? 
Even here, in a footnote to an appendix, one may not forget that the aim of these perturbative 
calculations is to understand, interpret, and improve the nonperturbative calculations. 
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Note that AA(0) = e'^^^l'^. 

With ali this machinery we can now state the main result of this appendix. 
The Feynman rules for vertices in standard references (e.g. ref. [ ^ for the 
Wilson action) are derived from the functional integrai, i.e. using '^[x). To 
obtain on-shell matrix elements, these rules must be supplemented by rules for 
contractions between iIjÌx) and conventionally normalized external states. They 
are 

• • • \\,^(t,x)---\---q{tp) • • •) ^ Ar(p)n„(e,p)e-^*+^^-, 



(•••gf(C,p)-- - I 



AA(p)n„(e,p)e+-^*-^J^, 



(C.23) 



(•••g(e,p) 



multiplied by the sign appropriate to the anti-commutation implied by the • • • . 
For ali states the momentum flow is physical, i.e. with (against) the charge flow 
for particles (anti-particles). 

For the specific action discussed in this paper, S'è, one finds (restoring a) 

^i(p) = C sin Pia 

(C.24) 

li{jp) = 1 + moa + ^rsCp^a^. 

The chromomagnetic and chromoelectric interactions do not modify these func- 
tions, but the kinetic corrections in Appendix ^ do. 

It is useful to record the small p expansion of the external line factor here: 



A^(p)^lat(^,p) = e 



-Afia/2 



2sinhMia 



8M2. 



u(C,0)+O(p^), (C.25) 



where the subscript "lat" abbreviates "lattice," and Mx is an "external line 
mass." For Sq 



1 



Mia? 



sinh^Mia + e^'-- ^^'^^^ 
For a unified treatment of fermions and anti-fermions in initial and final states, it 
is handy to note that eq. ( p.25D holds for positive and negative ^. The analogous 
expression for AA(p)?;iat(C,p) then follows from eq. ( |C.11| ). 

Unless rriQa <C 1 the lattice external line factor A/'(p)niat('^,p) deviates from 
the relativistic one. With our normalization conventions the relativistic analog 
of eq. (|a25|) is 

^ i-y-Tì 



^Urd(^,P) 



2m„ 



P 



n(e,0) +0(p^ 



(C.27) 



where the subscript "rei" abbreviates "relativistic." The rotations in sects. ^ 
and |A.3| are needed to convert the bracket of eq. (|C.25D info the bracket of 
eq. ( C.27| ), assuming niq = M2. 
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